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Abstract. The paper is concerned with mappings h: X y be- 
tween planar domains having least Dirichlct energy. The existence and 
uniqueness (up to a conformal change of variables in X ) of the energy- 
minimal mappings is established within the class (X, Y) of strong 
limits of homeomorphisms in the Sobolev space >^^'^(X, Y), a result 
of considerable interest in the mathematical models of Nonlinear Elas- 
ticity. The inner variation (of the independent variable in X) leads to 
the Hopf differential hzhzdz (8) dz and its trajectories. For a pair of 
doubly connected domains, in which X has finite conformal modulus, 
we establish the following principle: 

A mapping h € 3^2 (X, Y) is energy-minimal if and only if its 
Hopf-differential is analytic in X and real along dX. 
In general, the energy-minimal mappings may not be injectivc, in which 
case one observes the occurrence of cracks in X . Nevertheless, cracks 
are triggered only by the points in 9Y where Y fails to be convex. The 
general law of formation of cracks reads as follows: 

Cracks propagate along vertical trajectories of the Hopf differential 
from dX. toward the interior of X where they eventually terminate 
before making a crosscut. 
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1. Introduction 

Recently there have been new challenges and substantial work done [3, 4, 
5, 20, 28, 29] on minimizing the Dirichlet energy integral 

(1.1) ^^[h] = [[ \Dhf = 2 [[ {\h,f + \h,\^) dz, dz = dxi dx2 

subject to mappings h: X Y between bounded planar domains X and 
Y . The novelty of new problems is that the mappings in question are allowed 
to slip along certain parts of dH ; so-called traction free boundary problems. 
Motivated by mathematical models of nonlinear elasticity [7, 8, 9] the focus 
has been on finding a class of mappings as close to homcomorphisms as 
possible, in which the minimum energy is attained. The case when X and Y 
are simply connected is clear, by virtue of Riemann Mapping Theorem. Thus 
throughout this paper the domains X and Y will be predominantly multiply 
connected and of the same topological type, meaning that their complements 
C \ X and C \ Y consists of the same finite number of components, say 
i ^ 2. Prom these perspectives, the most desirable spaces in which to 
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look for the energy-minimal mappings are the class Y) C /^^'^(X, Y) 

of homeomorphisms and its strong closure J^(X, Y) in the Sobolev space 
/^^'^(X, Y), see Definition 2.4 for precise conditions. Note that the class 
^(X,Y) C y^^'^(X,Y) can be empty if dX contains isolated points as 
boundary components, but dY does not. Isolated points are removable 
singularities for homeomorphisms in y^^'^(X,Y). Such X can be viewed 
as a domain with punctures, which we call degenerate multiply connected 
domain. Apart from such a sitTiation we have: 

Theorem 1.1 (Existence). Consider a pair {X ,Y) of nondegenerate multi- 
ply connected domains, in which Y is a Lipschitz domain. Then there exists 
h:X^Y in the class jr2(X, Y) C W^''^lx,Y) such that 

(1.2) <^x[h] = inf{ ^x[f] : / e =^2(X, Y) } 

Specifically, we say that Y is a Lipschitz domain if every local arc in dY 
upon a suitable rotation becomes a graph of a Lipschitz function. 

Proof. We appeal to Theorem 1.1 and Remark 2.7 in [29] which tell us that 
the class ^jj^ (X, Y) of weak limits of homeomorphisms / : X Y in 
>ri'2(X,Y) (see Definition 2.4 ) and :^2(X,Y) C W^^^{X,Y) are equal. 
Armed with this result, the direct method in the Calculus of Variations 
concludes the proof. □ 

Theorem 1.1 calls for the following concept: 

Definition 1.2. Let X and Y be bounded domains. A map h G Jf2(X, Y) 
such that 

(1.3) <^x[h] = inf{ ^x[/] : / e ^2(X, Y) } 
will hereafter be referred to as an energy-minimal map. 

However the proof of uniqueness of the energy-minimal map is tricky and 
far more involved. 

Theorem 1.3 (Uniqueness). Let X be a nondegenerate doubly connected 
domain and Y a bounded doubly connected Lipschitz domain. Then the 
energy-minimal map h G J^2{X,Y) is unique up to a conformal change of 
variables in X. 

If one wants to find an easy and clear way of verifying whether a given 
map h : X Y is energy- minimal, one must look into the Hopf differential 
hzhz dz (8) dz . Here is the recipe. 

Theorem 1.4. LetX andY be bounded doubly connected domains, X being 
nondegenerate. Then a mapping h G Jif 2{X,Y) is energy-minimal if and 
only if its Hopf- differential is analytic and real along dX. 

It is surprising that this very useful criterion has not been established be- 
fore. For the convenience of the reader we devote the entire S3 to the inner 
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variations and the associated Hopf differentials; in particular, sec Defini- 
tion 3.8 that is pertinent to Theorem 1.4. The key observations in the proof 
of the uniqueness Theorem 1.3 is that the difference of two solutions to the 
same Hopf-Laplace equation has finite (not necessarily bounded) distortion 
on the set where at least one of the solutions has strictly positive Jacobian, 
see §9.1. 

Propagation of cracks is an interesting phenomenon not only in modern 
theories of elasticity and plasticity, materials science or microscopic crystal- 
lographic defects found in real materials but also from mathematical point 
of view. Let h G J^2i^,^) be an energy-minimal map between multiply 
connected domains. It should be noted that in general h, being a limit of 
homeomorphisms from X onto Y, has range in the closure of Y, /i : X — > Y. 
The fact, referred to as partial harmonicity, is that every h G J^2(X,Y)is 
a harmonic diffeomorphism of h^^{Y) C X onto Y, see [14] and Proposi- 
tion 6.4. There might, however, be sets in X which are taken into a single 
point in dY . 

Definition 1.5. Given a point a G dY , the term crack (or a-crack) in X 
refers to any connected component of the set { x G X: h{x) = a} . 

It will be seen (at least when Y is Lipschitz regular) via Lemma 2.2 that 
cracks always emanate from dX ; thus, are never reduced to a single point 
or any continuum in X . 

How to predict occurrence of cracks? One practical evidence is contained 
in the following result. 

Theorem 1.6. Let h G .y^2(X, Y) be an energy-minimal map, where X 
is a Jordan domain and Y a Lipschitz domain, both multiply connected. 
Suppose Y is convex at a boundary point a G dY , meaning that the set 
B CiY is convex for some ball B centered at a . Then a ^ /t(X) . 

This theorem may be interpreted as saying that cracks do not collapse 
into a boundary point at which Y is convex, see Theorem 12.2 and §13.2, 
see also [27] for similar phenomenon in higher dimensions. In spite of the 
impressive progress in the field, formation of cracks under energy-minimal 
deformations, to our knowledge, is not fully resolved. In this domain the 
best reference is the book by K. B. Broberg [12]. 

The problem clearly depends on the geometry of trajectories of the Hopf 
differential hzhg dz ^ dz. In this paper we give a detailed description of 
cracks in case of doubly connected domains. Accordingly, if the conformal 
modulus of Y (deformed configuration) is large relative to X then the Hopf 
differential hzhg dz dz is real and negative along dX. Consequently no 
cracks occur, even in the presence of points of concavity in the boundary of 
Y, see [20] and Theorem 7.1. If, on the other hand, the target is confor- 
mally very thin then the cracks are unavoidable. In this case the boundary 
components of X are horizontal trajectories along which the Hopf differen- 
tial is positive. The cracks are born in dX and propagate along the vertical 
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trajectories toward the interior of X where they eventually terminate. In 
particular, cracks create no crosscuts in X ; the domain X with cracks re- 
mains doubly connected, sec Theorem 7.3. 

A thoughtful analysis of the doubly connected case shows that one needs, 
in order to reveal the curiosities and the mechanics of cracks, investigate 
and explore the geometry of trajectories of the Hopf quadratic differential 
hzhg dz <^ dz by means of the level sets of h. Theoretical prediction of 
failure of bodies caused by cracks is a good motivation that should appeal 
to mathematical analysts and researchers in the engineering fields. 

2. Notation and Preliminaries 

2.1. Domains. We shall study a pair of bounded domains X , Y C ~ C 
of the same connectivity 1 ^ ^ < oo . This amounts to saying that each of 
the complements C \ X and C \ Y consists of £ disjoint closed connected 
sets. Let us introduce the notation, 

Xi, X2, X^ -the components of C \ X 

Yi , Y2 , . . . , Y£ -the components of C \ Y 
Their boundaries are exactly the components of dX and dY, respectively; 
8X1,8X2, dX( -the components of 8X 

(?Yi, 3Y2, dYi -the components of dY 

Note that adding any number of sets in {Xi,X2, ...,X^} to X results in a 
domain. 

2.1.1. Doubly connected domains. For a doubly connected domain X we re- 
serve special notation, X/ and Xo, for the bounded and unbounded compo- 
nents of C \ X , respectively. Thus 8X consists of two continua 8Xj and 
8X0 , referred to as inner and outer boundaries. 

2.1.2. Boundary Correspondence. Every homeomorphism h: X Y gives 
rise to a one-to-one correspondence between boundary components of X and 
boundary components of Y which, upon a suitable arrangement of indices, 
will be designated as 

(2.1) h:dX^^dY^, for = !,...,£ 

Precisely this means that the cluster set of dXi, under h is the bound- 
ary component dYi, . The components will be so ordered that the outer 
boundary of X corresponds to the outer boundary of Y . 
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2.2. The classes ^(X,Y) and J^d(X,Y). Let ^(X,Y) denote the class 

of all orientation preserving homeomorphisms h: X Y which preserve 
the order of boundary components in (2.1) . 

The following geometric concept of convergence [20] proves very useful. 

Definition 2.1. A sequence of continuous mappings /i j : X — >■ Y is said to 
converge cd-uniformJy to a mapping h: X — ;> Y if ; 

• hj ^ h , c-uniformly (uniformly on compact subsets of X ) 

• dist[hj{x), dY] dist[/i(x), dY] , uniformly in X 

cd 

We write it as hj — > h . 

Let J^rf(X, Y) denote the class of cd -limits of homeomorphisms in 
J^(X, Y) . Passing to the cd-limit one may lose injectivity but the bound- 
ary correspondence in (2.1) remains the same. It is a matter of topological 
routine to see that 

Lemma 2.2. For every h G ,^d(X,Y) it holds 

(i) Y C ^(X) C Y 

(ii) Given a boundary point a € 9Yjy C dY , every component of the set 
h~^{a) = { a; G X ; h{x) = a} (later interpreted as crack) touches 
the boundary component 5Xjy C 9X , i.e. its closure intersects . 

Another concept we are going to take from [20] is 

Definition 2.3. A mapping ^: X — >■ Y is called deformation if 

• /iG>ri'2(x) 

• J{x, h) ^ almost everywhere 
. JJ^Jix,h)dx^\Y\ 

• There exist sense preserving homeomorphisms fj : X Y (called 

an approximating sequence) converging cd -uniformly to h . 

Let us emphasize that the approximating homeomorphisms fj need not 
have finite energy. The set of deformations is denoted by 33(X, Y) . 

2.3. The Dirichlet Energy. The Dirichlet energy (or Dirichlet integral) 
for functions h:X^C in the Sobolev space /^-"^'^(XjC) is defined by 

(2.2) ^^[h] = [[ \Dh\^ = 2 [[ (|/i,|2 + \h,\^) dz, dz = dxi dx2 

Here and throughout, we take advantage of the complex variable z = xi+i X2 
and partial derivatives (Wirtinger operators) 

^ _dh _\f dh^ _.dh\ ^ _dh _ \ f dh ^ -dhX 
^ dz 2 \dxi 8x2 J ' ^ dz 2 \dxi 8x2 J 
Thus, the Jacobian determinant is 

(2.3) J{z, h) = Jh{z) = det Dh{z) = \h^f - \h^f 
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2.4. The Space J^(X,Y) and other Relevant Classes of Mappings. 
Various classes of mappings of finite energy will be considered. 

Definition 2.4. Here they are: 

• j^{x, Y) = jf{x, Y) n /r^'2(x, Y) . 

• ^(X, Y) is the closure of (X, Y) in strong topology of W^''^ (X, Y) 

• ^(X,Y) is the closure of ^(X,Y) in weak topology of W^'^{X,Y) 

• Diff^ (X,Y) consists of <^°°-diffeomorphisms in ^(X,Y) . The fact 
is that :^(X, Y) = Diff^ (X , Y) , see [23]. 

• J^}^{X,Y) stands for the family of all weak limits of homeomor- 
phisms hj:X°-^Y , j = 1,2, ... , in the Sobolev Space W^''^{X, Y). 

Let us call attention to an important detail that h , in any of the above 
classes, has a range in the closure of the target domain; that is, ^: X ^ Y . 

Remark 2.5. In general, the family of all weak limits of a set in a Banach 
space need not be weakly closed. That is why we have only the inclusion 

^^^(X,Y) ^ J^(X,Y). However, if X is a Jordan domain and Y is a 
Lipschitz domain, both not simply connected, then 



(2.4) ^i;;'(x, Y) = jr,(x,Y) = jr2(x, Y) = ^2(x, y) c =^cd(x, y). 



For the first two equalities we refer to [29]. The symbol Jf2(X,Y) de- 
notes the class of strong limits of homeomorphisms up to the boundaries, 
hji X Y , in the Sobolev space ^'^'^(X, Y) ; the equality is immediate 
from [24]. The last inclusion is then obvious. 

Finally, we infer from [20] that 

Lemma 2.6. Let X and Y be multiply connected domains, X being non- 
degenerate. If a sequence {hj} C 2)(X, Y) converges weakly in W^'^{XY) , 
then its limit belongs to S)(X, Y) . Thus, in particular 



2.5. Monotone mappings. Monotone mappings were born in the work of 
C.B. Morrey [42]. We will be concerned with compact subsets of so let 
us combine Morrey's original idea with a theorem of G.T. Whyburn , see 
page 2 in [41], to make the following definition: 

Definition 2.7 (monotonicity) . A continuous mapping /: X Y be- 
tween compact metric spaces is monotone if for each connected set C C Y 
its preimage /~^(C) is connected in X. 

Theorem 2.8 (Kuratowski-Lacher) . Let X and Y be compact Hausdorff 

spaces, Y being locally connected. Suppose we are given a sequence of mono- 
tone mappings fk ■ X Y converging uniformly to a mapping / : X — )• 
Y , then f : X "^^^ Y is monotone. 

See [36, 37] and [47, 52] for further reading about monotone mappings. 



(2.6) 



(2.5) 



J^i;;'(X,Y) C 2)(X,Y) C jrcd(X,Y) 

Y C h{X) C Y , for every h £ ^^^(X, Y) 
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2.6. Equicontinuity and monotonicity. 

Lemma 2.9. Let X and Y be Lipschitz domains of connectivity i ^ 
2. Then, every h G J^l^{X,Y) extends to a continuous monotone map 
h : X Y . The boundary map h : dX dY is also monotone. More- 
over, 

(2.7) \h{xi) - hix2)f ^ ^^^^^^—^[[\Dh\\ xi,X2GX. 



log(e + ^) 



Proof. It is known that a homeomorphism h: X Y between Lipschitz 
domains in the Sobolev space 'W^''^{X,Y) extends continuously up to the 
boundaries [28]. And it is topologically clear that a continuous extension of 
a homeomorphism h : X Y results in monotone mappings h : X Y 
and h : dX dY . These properties carry over to the uniform limits of 
homeomorphisms. Thus it only remains to justify the uniform estimates 
(2.7). Local estimates like this are well known for monotone mappings of 
Sobolev class ■W^'^(X,Y) [21, 26]. The global inequality follows by the 
standard method of extending h beyond the boundaries. It is at this point 
that the Lipschitz regularity of the domains is required again. We leave the 
routine details to the readers. □ 

Remark 2.10. Similar arguments, including the extension procedure, apply 
to a pair of simply connected Lipschitz domains if the mappings are fixed at 
some interior point ( h{xo) = yo for some Xo G X and yo G Y) or at three 
boundary points ( h{xi) = yi , h{x2) = y2 , h{x2) = y2 )• We shall return to 
this later case in §12.1, where we will discuss quadrilateral mappings. 

2.7. Rado-Kneser-Choquet theorem. The following proposition strength- 
ens the well know theorem of Rado-Kneser-Choquet. 

Theorem 2.11. A harmonic map h: Q ^ C of a Jordan domain, which 
extends continuously as a monotone map of dO, onto a boundary of a convex 
region is a -diffeomorphism of J7 onto this region. 

For the proof and the best general references see [15]. 



3. Inner Variation and Hope Dieferentials 

Although the concept of inner variation has been used in the Calculus 

of Variations for a long time [11, 13, 32, 50, 49], see also more recent ones 
[14, 25, 43, 53], some of its nuances are still to be scrutinized. 

3.1. Inner Vziriation. Throughout this section X cM? ^ C is a bounded 

domain whose points, also called variables, will be denoted by x = xi + 1x2 ■ 
To cover all the instances discussed later, we distinguish a compact subset 
IK C dX , possibly empty or the entire boundary. 
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3.1.1. Change of variables. The term change of variables in X that is fixed 
at the compact subset IK C dX , refers to a "^""-diffeomorphism : X X 
which is continuous up to dX and satisfies ^{x) = x for x e K. 

3.1.2. Variation of variables in X . This is a one-parameter family of change 

of variables {^''^}„£^<£<£^ , such that 

(i) The function {e,x) "^'^{x) is ^°°-smooth in {—So , £o) x X 

(ii) = id in X 

(iii) = id in K, for every e G (— £o , So) 
Define A = and note that DX=^\ ^D^^^ . 

de I £=0 de I £=0 

3.1.3. Inner variation of a mapping. Let h he a function in the Sobolev 
space ^^'^(X, C) and {^''} a variation of variables in X that is fixed 
at K. C dX . The family h^ = h o is called an inner variation of 
h . Observe that all the mappings h^ have the same range as h , which 
is one of the desired properties that motivates the use of inner variations. 
Additional assumption on the behavior of ^'^ near dX will be needed in 
order to ensure that G >^^'^(X, C) . In our appHcations, however, this 
property will always be satisfied, so there is no need for formulating any 
explicit conditions here. Nevertheless, the interested reader may notice that 
if all the diffeomorphisms \1' ^ : X X are X-quasiconformal; that is, 
\D^^\ < KdeiD^^ , then 



// \Dh^{x)\^dx [[ \Dh{x) 



' dx < oo 



^ def d$^ 



de 



£ = 



3.1.4. Critical points. Choose and fix a variation {^''^}_£^<£<£j^ of variables 
in X. We say that h G #'^'^(X, C) is a critical point for {^^} if 

de £=0 

Lemma 3.1 (Integral form of the variational equation). Every critical point 
for {'^^}-eo<e<eo Satisfies the equation 

(3.1) Re J J hzhzXzdz = 0, where 

Proof. The peculiarity of our derivation of (3.1) lies in pointing out the 
relevance of the equation to the Cauchy- Green stress tensor of h and the 
AhlfoTS infinitesimal deformation operator for A . These connections were 
certainly overlooked in the literature; both are worth noting for prospective 
generalizations. For example, the same ideas work for the conformally in- 
variant 77,-harmonic integrals in M" , see [27] . 

Let us begin with the composition rule Dh^ = Dh{^^) ■ . We express 
the Hilbert-Schmidt norm of the differential matrix Dh^ via scalar product 
of matrices 

= {Dh^ I Dh') = {Dh{^^) ■ D*^ I Dh{^^) ■ D^^) 
^^■^^ = ( D*h{^^) ■ Dh{^^) I D^^ ■ D*^^ ) 
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where D* stands for the transposed differential. Then we integrate it over 
X by substitution ^ = 



(3.3) <fx[^1 = jj (^D*h-Dh 



dct 



At this point one may recall the Cauchy-Green stress tensor of h 

C[h] = D*h-Dh 
as well as Weyl's conformal tensor of 



det D^^ 



Since = id , 



d*i| = A and 



= DX , the infinitesimal varia- 



tion of W[^'^] at the identity map is readily computed as 
dW[*^ 



de 



£=0 



DX + D*X - (TV DA)/ 



Here too, wc recognize the Ahlfors infinitesimal deformation operator [1, 2, 
48] which defines the conformal component of the matrix DA , 



SA = ^ [DA + D*A - (TV DA)/] 



Re Xz , Im A^ 
Im A2 , — Re Xz 



The anticonformal component is afforded by the Ahlfors ' adjoint operator 



AA = ^ [DA - D*A + (TV DA)/] 



ReA^ , 
Im Xz , 



— Im A^; 
ReA^ 



Thus we have a pointwise orthogonal decomposition DA = SA + AA . Now, 
a straightforward computation gives the variation of the energy integral 



(3.4) 



d^x[^ 
de 



£=0 



D*h ■ Dh 



SA)dC 



= 8 Re / / hzhz Xz dz 



because {D*h ■ Dh | SA) = 4Re{hzhzXz) 



□ 



3.2. Hopf-Laplace equation. Every complex function A G '^q°°(X, C) 
gives rise to a variation '^^{z) = z + eX{z) of variables in X that is 
fixed at every point of 9X . This variation is legitimate for any function 

h G /^^'^(X, C) that is fixed at a given compact subset K C 9X . In par- 
ticular, applying (3.1) to such A we infer via the classical Weyl's lemma 
that 
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Proposition 3.2 (Hopf-Laplace Equation). A function h G C) 
that is critical for all inner variations (subject or not to boundary con- 
straints ) must satisfy the Hopf-Laplace equation in X , 

(3.5) -qZ hz hz = , (in the sense of distributions) 

equivalently, 

(3.6) hzhz=(p{z) , for some analytic function <^€^^(X) 

If one allows the critical mapping h to slip along an arc 7 C 9X then an 
additional equation on 7 will emerge, which in turn yields a specific form 
of the analytic function . Let us begin with an example of such situation. 

Example 3.3 (Critical solutions in an annulus). Consider a traction free 
problem in an annulus X = {z; r < \z\ < R}; that is, allow h: dX dY 
to slide along the boundary circles of X . 

Proposition 3.4. The Hopf-Laplace equation for a function h G #'^'^(X,C) 
that is critical for all inner variations in X takes the form 

(3.7) hzkz=-^, for all zeX 
where c is a real number 

Proof. We expand the Hopf product of h into a Laurent series 

00 

hz hz = (p{z) = ^ GnZ" 

n=—OD 

Then wc test the integral equation (3.1) with the following variations of 

variables in X , 

\l — e'^{az'^ — az'^Y]^!'^ 

where a can be any complex number and e any sufficiently small real num- 
ber. Evidently |*^(2;)| = \z\ , so X X can easily be shown to be a 
-diffeomorphism. A short computation gives 



A 



def d^-^ 



= z[az — az 

e=0 



k „^k\ X= = kazz'^ ^ 



de 

Put these values of Ag into the equation (3.1) to obtain 

= Re // ( V On 2") kazz^-^ dz = Re{kaak-2) II kP^"^ d^ 

for every complex number a . This yields afc_2 = , except for fc = . We 
just proved that ^p(z) = a_2 z~'^ . To see that the coefficient a_2 is real 
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we test the integral equation (3.1 again, but with the following variation of 
variables in X . 

^ ^ [l + e2iog2|^|]i/2 de 
Then we find that 



= izloglzl and = i- 
e=o z 



= Re // i- dz = — Im(a ,) // \z\ ^ dz , hence a 

JJx z yjx 



2 e 



as desired. □ 

3.3. Variations along a boundary arc. In the study of traction free elas- 
tic deformations h: X Y one prescribes the values of /i in a compact 
set K C 5X while allowing to slide along the rest of the boundary. The 
critical points, subjected to mappings which are traction free along the set 
r = 5X \ K , must satisfy not only the Hopf-Laplace equation (3.6) in X 
but also an equation on F . We begin with a smooth Jordan arc in V , and 
later extend this concept to more general settings, sec Definition 3.8. 

3.3.1. Boundary arcs. A subset 7 C 9X is said to be -smooth bound- 
ary arc if there is an open set U C C and a -smooth diffeomorphism 
F:U°-^C such that 

• U n 5X = 7 , F(7) = M 

• F{Xnl!) = C+ = {z = xi + ix2; X2 > 0} 

We say that F : U C transforms 7 into a line segment in U , and refer 

to U as vicinity of 7 . 

3.3.2. Test functions. A complex function A € '^o^(U) is a test function fox 
the arc 7 C 9X if the linear differential X{z) dz is real along 7. This means 
that for any (or just one) parametrization of 7 , say 7 = {z{t) : a <t < 
b }, it holds \{z{t)) z{t) G M. Here, as usual, z{t) denotes the derivative 
of z{t) . In other words, the complex number X{z) represents the tangent 
vector at z G 7 . A diffeomorphism F : U C which transforms the arc 
7 into a line segment in M gives a function A o F~^ G ^o'(C) that is real 
along M. 

Lemma 3.5 (Existence of the variation). Given a smooth boundary arc 
7 C 9X and its test function X G Co(U) , there exists a variation of variables 
*=:X^X along -f , such that ^I^^q"^^ =A 

In particular, by Lemma 3.1, we have 

Corollary 3.6. If h & y^^'^(X, C) is a critical point for all variations along 
7 then, 



(3.8) ^^JJ ^^^^^^^^ = 

for all test functions A G Cq (U) whose linear differential A dz is real on 7 . 
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Proof. An explicit construction is as follows. Fix a -diffeomorphism 
F : U C which transforms 7 into a line segment in M . Then define the 
desired one-parameter family of diffeomorphisms by the rule, 

*^ = (F + eXF^ + eXF^) : U n X ^ C 

with |e| < £0 sufficiently small, and set ^^(z) = z in X \ U. 

Here Fz and denote the complex derivatives of F : U C . □ 

With Corollary 3.6 in mind, we introduce the following concept 

Definition 3.7. A holomorphic quadratic differential (p{z) dz (Si dz, with 
if € ^^(X) , is real along a smooth arc 7 C dX , in a weak sense, if 

(3.9) Re JJ <p{z) Xz dz = 

for all test functions A G Cg (U) whose linear differential A dz is real on 7 . 

This definition is not more general than the classical one, but is convenient 
when dealing with critical points of inner variations along 7 . Indeed, it turns 
out that ip admits continuous extension up to X U 7 and as such fulfils the 
classical condition: 

(3.10) cp{z{t))iitf gM 

for some (also for all) smooth parametrization of 7 , ■y = {z{t) : a < t < b } . 
In full generality one might say that: 

Definition 3.8. A holomorphic quadratic differential (^(x) dz dz , with 
(p G =5f ^ (X) , is said to be real along dX if upon conformal transformation 
of X onto a "^^-smooth domain the differential becomes continuous and real 
along its boundary. 

The definition is legitimate, because the choice of conformal transforma- 
tion plays no role. 

3.4. Quadrilateral mappings. In quasiconformal geometry the term quadri- 
lateral is used to describe a simply connected Jordan domain Q C C to- 
gether with the ordered set of four points a,b,c,d G listed counterclock- 
wise along the boundary and called corners of Q . The corners disconnect 
dQ into four open Jordan arcs; two of them, having no common endpoint 
are called horizontal sides and the other two vertical sides. The exemplary 
quadrilaterals, from which other quadrilaterals are derived via conformal 
transformations, are the rectangles. But we shall choose annular sectors as 
representatives of quadrilaterals. 
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3.4.1. Annular Sectors. Given ^ a < /3 ^ 27r and < r < R , the 
annular sector 

^ = ^^{r, R) = {z;r <\z\<R, a<argz</3} 

is viewed as a quadrilateral whose corners listed counterclockwise are a = 

re* °' ,h = i?e* ° , c = Re^ ^ , d = re* ^ . The circular arcs in 5 S are viewed 
as vertical sides and the rays as horizontal sides. 



3.4.2. Quadrilateral Maps. 

Definition 3.9. A homeomorphism h: X Y between two quadrilaterals 
which has well defined limits at the corners of X and that these limits 
coincide with the corners of Y (in the respective order) is called quadrilateral 
map. 

The reader is cautioned that the inverse homeomorphism : Y X 
may not be a quadrilateral map. It should be noted, however, that the clus- 
ter set of any side of X under a quadrilateral map h : X Y is exactly the 
corresponding side of Y . We have a quadrilateral conformal transformation 
z j-^ia+p z-az q£ rectangle TZ = (0, L) x (0, 1) onto the annular sector 
S = S^(r, R), L= ^logf . 

Minimization of the Dirichlet energy among quadrilateral mappings h : § 
Q calls for inner variations : S S that keep the corners of the sector 
fixed, while allowing slipping along the sides. Call them quadrilateral varia- 
tions. With a few additional tricks the same inner variations as are used in 
Proposition 3.10 show that 

Proposition 3.10. The Hopf-Laplace equation for a quadrilateral mapping 
h G >^^'^(S,C) that is critical for all quadrilateral variations in S takes the 
form 

(3.11) (f>{z) = h^k^ = ^ , for ah z G S 

where c is a real number 

TT — Q 

Proof. First observe that the power conformal map x{^) = c*" z takes 
the sector 'E>a{r, R) onto an upper half annulus Sg ('^'j -^0 ' ^'^d it is a quadri- 
lateral map. Another useful observation is that the Hopf-Laplace equation 
(3.11) is invariant under power type transformations of sectors. It is there- 
fore enough to proof (3.11) for an upper half annulus. Quadrilateral varia- 
tions that we are going to explore this time are: 

so 

A 



de 



= zz'^ , \z =kzz^ ^ 

£=0 
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Curiously, the counterpart of = is 



z 



T£/ X 1 + ^Arg 2; ^ d$' , , . 

[1 + e^kig'^zYi^ de £=0 z 

though we can dispense with this variation. But we need a Laurent series 
expansion 

00 

(3.12) hzhz=(p{z)= ^ a„z" , an = an 

n=—oo 

The natural way to obtain such formula is by extending ip to the whole 
annulus A = S+ U S_ , where S+ = S and S_ is a reflection of § to the 
lower half plane. Wc extend cf) by setting (f){z) = (l>(z) for Im^ < 0. It 
is immaterial how (j) is defined for Im z = . Clearly, the extended func- 
tion, still denoted by (f) , lies in =Sf^(A) . Precisely, ||0||^i(a) = 2 ||^||^i(s) ■ 
Let r] G ^^(A) be any test function, to be used for verifying that the 
distribTitional Cauchy-Rieman equations in A are satisfied. We can write: 

Re / / (j){z) % dz = Re / / 0(z) r]z dz + Ke (j){z) rfz dz 
JJa J Js+ J Js- 

= Re / / (j){z) r]z{z) dz + Re (l){z) r}z{z) dz 

= Re / / (p{z) riziz) dz + Re 4>{z) riz{z) dz 
J J§+ J Js+ 

= Re ff (j){z) [r}z{z) + ^zi^ ] dz 

J Js+ 

= Re (l){z) \z{z) dz = Q 



where \{z) = ri{z) + tjiz) belongs to and is obviously real on M. 

The last equality follows from Lemma 3.5 . Replacing 77 by z 77 we conclude 
that J j'^ (piz) rjz dz = , which means that cf) is analytic in S , by Wcyl's 
Lemma. Having the Laurent expansion at hand we now test the integral 
form of the variational equation (3.1) , 

Re // hzhzXzdz = 0, with =kzz''~^ , A; = ±1, ±2, ... 



lis 



To this effect note that Re / Jg a„ 2;" dz = a„ Re / /g z"' Xz dz = , 
except for the case n = k — 2. We just proved that (j){z) = a_^ z~'^ , where 
o_2 = a_2 . That a_2 must be real can also be seen by testing the variational 



equation with A^ = z ^ 



z 

□ 



3.5. Vertical and horizontal arcs of a quadratic diflferential. Let 

(^(z) dz ® dz be a holomorphic quadratic differential in X . A vertical arc 
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is a '^°°-smooth curve 7 = 7(t), a < t < b, along which 
[7(t)]V(7(i)) < 0, a<t<b. 

A vertical trajectory of in X is a maximal vertical arc; that is, not properly 
contained in any other vertical arc. In exactly similar way are defined the 
horizontal arcs and horizontal trajectories, via the opposite inequality. 

If X is a circular annulus A = A{r, R) and ip{z) dz <S> dz is real along 
its entire boundary then (p{z) = cz^'^ , for some c G M . For c > the 
concentric circles Cp = {pe'^ : ^ 9 < 27r} , p E [r,R] , are the vertical 
trajectories, whereas the rays TZ^ = {pe^^ : r < p < R} , 6 € [0, 27r) , are 
the horizontal trajectories. For a negative c, this holds in reverse order. 
These two cases exhibit different behavior in regard to the formation of 
cracks. 

Every continuous mapping h: X — > has well defined multiplicity func- 
tion which is measurable so one can speak of the essential supremum. We 
are concerned with mappings of bounded multiplicity, meaning that there 
is 1 ^ M < 00 so that 

(3.13) #{ X G X ; h{x) = y} ^ M , for almost every y G M^. 

Indeed, for h G J^f^(X,Y) we have 

rc G X ; h(x) = y} ^ 1 , for almost every y G M'^ 
see [20, Lemma 3.8]. The following Proposition is from [14, Proposition 5.1]. 

Proposition 3.11. Suppose h G '^(X, C) Pi Will}{X,C) satisfies the condi- 
tion (3.13) and solves the Hopf-Laplace equation 

hzhz = f , where f ^ is a holomorphic function in X C C 

Then for each yo G the union of all vertical trajectories of the quadratic 
differential ip{z) dz ® dz that intersect h^^{yo) has zero measure. 

3.6. Bizarre solutions. Complex harmonic functions in the Sobolev space 
>^^'^(X,C) are among the solutions of the Hopf-Laplace equation. But 
there are many more, sometimes surreal ones. Harmonicity is lost exactly 
at the points where h fails to be injective [22]. Consider, for example, an 
origami folding of a square sheet of paper, denoted by X . Upon a countable 
number of suitable folds all the points in dX stock up into one point, say 
the origin. In this way one obtains a piece-wise orthogonal Lipschitz map 
h : X Y which vanishes on 9X . In fact the range of the differential 
Dh may consist of a finite number (eight is enough) of orthogonal matrices. 
Therefore, at almost every point z G X we either have hz = or hg = . 
The interested reader is referred to [30, 34, 35] for more details. Thus we 
have a nonzero solution of the homogeneous Dirichlet problem 

r hzhi = , h£Lip(X, Y) 
\ /i = , on ax 
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There are, perhaps, even more exotic solutions to the homogeneous Hopf- 
Laplace equation, hzhz = , interesting not only in their own right. 

3.7. The natural domain of definition and admissible solutions. It 

becomes clear that in order to build a viable theory of the equation (3.1) 
one must assume that the Jacobian determinant is nonnegative almost ev- 
erywhere. Such reflections urge the following concept. 

Definition 3.12. The class #'+'^(X,C) of mappings in W^l'^{X,C) whose 
Jacobian J{x, h) ^ almost everywhere in X will be regarded as the natural 

1 2 

domain of definition of the Hopf-Laplace operator. A map h G W+'^{X, C) 
which satisfies the equation J| [hzhz) = in the sense of distributions, will 
be referred to as admissible solution. 

Critical points for the Dirichlet energy within any of the classes in Defini- 
tion 2.4 comply with these requirements. It should be remarked that all the 
admissible solutions to (3.1) are locally Lipschitz continuous [14] and [25], 
but not ^^-smooth in general. See also Theorem 12.3 which, in a specific 
case, asserts Lipschitz regularity up to the boundary. 

4. Harmonic Replacement 

The main idea in the proof of Theorem 1.6 is to show that every h G 
J^^^{\,Y) = M^2{^i^) can be modified in a neighborhood of the set 
h~^{a) C X, so-called a boundary cell, with lower energy, unless h al- 
ready takes the boundary cell diffeomorphically into Y . To this effect we 
need a thorough analysis of harmonic replacement procedure. 

Let h e J^^l^(K,Y) and 5 C C be a ball. Assume that the set D = 
B n Y is nonempty and convex. Then, because of monotonicity, the set 
U = {a; G X ; h{x) G i? fl Y} , is an open simply connected domain, called 
a cell in X . More specifically, we call U an inner cell if i? C Y and a 
boundary cell otherwise. 

Proposition 4.1. LetY be a Lipschitz domain and h G J^j|^^(X,Y) . Then 
to every cell U C X there corresponds (unique) mapping G Jf'2(X, Y) 
such that 

(i) h^ = h:X\U°-^Y\B 
(a) /ip : U D is a harmonic diffeomorphism. 
(tti) ^x[/i„] ^ ^x[h] 

Equality in (Hi) occurs if and only if h = h^. We refer to as harmonic 
replacement of h in U . 

Before passing to the proof let us make a few useful observations. First, 
Proposition 4.1 is not affected by a conformal change of variables within the 
domain X . Thus it involves no loss of generality in assuming that X is 
a Schottky domain (bounded by circles), so Lipschitz regular. Under this 
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assumption, equalities in (2.4) apply. In particular, we dispose of a sequence 
of homeomorphisms hj : X Y which converge to h uniformly on X and 
strongly in ^i'2(X,Y). The limit map h extends continuously up to the 
closures, which we continue to denote by /i : X Y . The extended map 
h : X Y and the boundary map h : OX dY are monotone. 

Proof. The case of an inner cell is particularly convenient to begin with. 
Later, with additional topological arguments, we will be able to exploit this 
proof for the boundary cells as well. 

Case 1. The inner cell. Suppose B C Y so U = h-^{B) C X. bmce 
h is monotone, it follows that U is a simply connected subdomain of X. 
Consider an increasing sequence Bi C B2 C . . . , UB„ = B, of balls and the 
corresponding cells Un = h~^{Bn) 

(4.1) U„_i c U„_i c U„ c U„ c U„+i 
and note that 

dUnCh-HdBn) = rn. 

Here r„ is a continuum disconnecting C into two components; the bounded 
component equals U^. Thus U^-i is a compact subset of the inner comple- 
ment of r„ . We look at the induced cells Un = hJ^{Bn) , j = 1,2, ... . 

Their boundaries = dVi, = hJ^{dDn) are closed Jordan curves. Since 
hj converge uniformly to ^ , it follows that 

hm sup {dist(x, r„) ; x G F^} = 

j-5>00 

which in turn implies that for sufficiently large j, say j ^ j„, we have 
C Un+i and U„_i is contained in Un -the bounded component of C \ 
r^. We now appeal to Rado-Kneser-Choquet theorem, see Theorem 2.11. 
Accordingly, we may extend the boundary homeomorphism hj^ : Fn" — t- dBn 
inside the cell Un" to obtain a homeomorphism hj^ : X Y such that 

. hj^ = hj,^:X\V(r-^Y\Bn 

n is a harmonic difFeomorphism 

• <%[^i„] ^ '%[/ij„] =^ M , this is a bound independent of n. 
Passing to a subsequence of {hj^}, if necessary, we define h G J^j^''(X,Y) 
to be a weak limit of {hj^}. This sequence also converges uniformly to 
h: X Y , because of the uniform bounds in (2.7). Moreover, 

h = h: X\U°-^Y\B . 

Let us demonstrate that 

(4.2) h: U 2^ B is a harmonic diffeomorphism. 

To this effect, fix any sequence {KnjJ^i of subdomains compactly con- 
tained in Un , such that UKn = U. Since /ij^ : Un — >■ C, k ^ n, converge 
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uniformly to h: K.^ — )• C we shall infer that either h: ¥^ — )• C is a diffeo- 
morphism or J(x, /i) = in at least for large values of n. This fact 
is none other than a harmonic variant of the classical Hurwitz theorem for 
analytic functions. 

Lemma 4.2 (Hurwitz Lemma). If a sequence of harmonic homeomorphisms 
fk - ^ ^ C converges c-uniformly to f : ^ C then either f is a har- 
monic homeomorphisms (actually '^'^ -diffeomorphism) or its Jacobian de- 
terminant vanishes identilcally, J{x, /) = . 

The proof of this lemma presents no difficulty; simply, consider a sequence 

{^^} of analytic functions. Curiously, it also holds for p-harmonic map- 
pings, but in this case the proof requires much more work, see [29]. 

Returning to our demonstration, the equality J{x, h) = in is ruled 
out by the following computation 

// J{x,h)dx= lim // J{x,hji.) dx ^ lim // J{x,hj^)dx 
J Jv ^-^'^ J Jv J Jvi'' 

k 

= lim \hj,{Vf)\ = lim |5,J = \B\ > . 

Here, in the first equation, we appealed to the =Sfj^^, - weak convergence of 
nonnegative Jacobians, a property discovered by S. Miiller [44], also see [26, 
Theorem 8.4.2]. Thus, for large n we have JJ^ J{x,h)dx > and so 

J{x,h) ^ in IKjj . M^e infer that /t is a local diffeomorphism in U. Since 
h is a c-uniform limit of homeomorphisms, the map h : U B must be 
a global homeomorphism. Finally, we apply the Dirichlet principle which 
asserts that if in a bounded simply connected domain U two functions 
h,he -^(U) n W^'^iV) coincide on dV , then 

'%[^] ^ ^v[h] ) whenever h is harmonic 
Equality occurs if and only ii h = h onU. The proof of Case 1 is complete. 

Remark 4.3. In the above proof wc actually did not use any regularity of X 
or Y. Thus Proposition 4.1 holds for arbitrary domains X and Y , provided 
B cY ; that is , U = f~^{B) is an inner cell. It is the case of the boundary 
cells that we really benefit from Y being a Lipschitz domain. 

Case 2. The boundary cell. Apart from a few adjustments, necessi- 
tated by some geometric complications near the boundary of X , the proof 
differs very little from that in Case 1. The idea is to extend the mapping 
h: X Y beyond X an Y so the boundary cell of h will be enlarged 
accordingly to become an inner cell of the extended mapping. In all that 
follows X will be a Schottky domain (i.e. bounded by circles) and Y a 
Lipschitz domain. 
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Lemma 4.4 (Extension Lemma). There exist neighborhoods D X and 
Y+ D Y and an extension operator 

so that for every h € J^2(X,Y) its extension h G Jf'2(X+,Y+) satisfies : 

(i) h{x) = h{x) , for x G X 

(ii) /i: X+\X^^ Y+\Y 
(iii) 

// \Dh{x)fdx ^ C„ f f \Dh{x)\^ dx 
J Jx+ J Jx 

(iv) Let a sequence of homeomorphisms hj : X Y converge to h 
uniformly and strongly in W^''^(K,Y). Then the extended map- 
pings hj : X+ Y-|_ are also homeomorphisms; they converge to 
h uniformly and strongly in >^^'^(X+ , Y_|_) . 

Proof. A brief sketch of such extension is the following. First, with the aid of 
a bi-Lipschitz automorphism of C we transform the target domain Y onto 
a circular domain, say Y' . There is no need for such a transformation of X 
because it is already a circular domain, record this assumption as X' = X . 
Next we reflect X' and Y' about their boundary circles and add the reflected 
images to X' and Y' , respectively. We obtain the neighborhoods X^ D X' 
and Y^ D Y . Then we return, via inverse bi-Lipschitz automorphism, to 
the desired neighborhoods X+ D X and Y+ D Y . These same reflections 
give rise to the extensions of the mappings hj , denoted by hj : X+ 
Y_|_ , and of the limit map h , denoted by h : X-|_ Y_|_ . The interested 
reader is referred to [29] for details, especially to rather unexpected nuance 
about continuity of the composition of a Sobolev mapping with bi-Lipschitz 
transformation of the target domain. □ 

Armed with Lemma (4.4) we are now able to repeat most of the arguments 
in Case 1. Recall the boundary cell U = {x € X; h{x) € n Y} for 
h . Here we may, and do, assume that the ball B lies in X-|- , so the set 
{x G X; h(x) G } is an inner cell for h . For convenience, we continue 
to use the same notation for the induced cells as in Case 1; namely, 1[J„ = 

h-^{Bn) , = hj~^{Bn) and Fi = dUi . These are subsets of X+ . Fix 
any sequence {IK„} of domains compactly contained in H X , such that 
IJ^2 = UnX. However, this time harmonic extensions of the boundary 
homeomorphisms hj : dBn are useless. They do not take the arcs 

UiiCidX onto the corresponding arc S„n5Y ; even if they do, such harmonic 
extensions may not coincide with h on these arcs. That is why we should 
perform harmonic replacements only on a portion of U^^ ; precisely, on the 
set Un n X . It is important to observe that n X is a simply connected 
Jordan domain. To see this, consider the open Jordan arc 7 = S„ fl dY , 
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and its homeomorphic preimage P = hj^{-f) cUindX. The endpoints of 

the arc /3 he in the boundary of Un , as is easy to demonstrate. Thus /3 is a 
crosscut of the simply connected domain Un • This is surely an elementary 
topological fact that such a crosscut /3 splits Un into two simply connected 
subdomains. The subdomain of interest to us is exactly UnflX . Now, upon 
harmonic replacements of hj in Un n X, the proof of Proposition 4.1 runs 
as in Case 1 with hardly any changes. □ 

5. Proof of Theorem 1.6 

The case when X and Y are simply connected is obvious; in this case h 
is a conformal mapping. Thus let X and Y be multiply connected Jordan 
domains, Y being a Lipschitz domain. Consider the boundary cell U = 
h~^{Br\Y), where a G BnY , and the harmonic replacement hu € c^2(X, Y) 
as in Proposition 4.1. Since h is energy-minimal, we have £x[h] ^ iSx[^u] • 
Therefore, by condition (iii) we must have equality £x[h] = ^x[^u] which 
yields that h = h\}. The proof concludes by observing that 

h:X\V ^Y\B , thus a /i(X \ U) 

h:V°-^Ii, thusa0/i(U) 

Therefore, a /i(X \ U) U /i(U) = /i(X) , as claimed. □ 

6. The CLASS J^^^(X,Y) revisited 

We emphasize that the boundary points of Y may be in the range of 
h G J^^;(X, Y) , which is the major source of difficulties. Nevertheless, by 
Lemma 2.6 we always have 

(6.1) Y C ^(X) C Y, whenever h G J^^f (X, Y) . 

1 2 

Lemma 6.1. For h G ^^ji^ (X, Y) the Jacobian determinant vanishes al- 
most everywhere in X \ h~^{Y) = h~^{dY) = {x: h{x) G dY} . 

Proof of Lemma 6.1. Consider the set B = X\ h-^{Y) = h'^idY) and let 
K be any compact subset of B. We invoke the sequence of homeomorphisms 
hj : X Y converging to h weakly in the Sobolev space 1^^'^(X, Y). By 
■^loc ■ weak convergence of nonnegative Jacobians, [26, Theorem 8.4.2] , we 
can write 

/ J{x,h)dx= lim / J{x,hj){x) dx = \im.\hj(K)\. 
Jk Jk 

Claim. We have hj {K) C Y^ , where 

Yj = {y eY: dist(2/, dY) ^ ej} and ej = sup||/fe - /||^oo(k) ^ 0. 
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Indeed, let y = hj{x) G hj(K) for some a; € K, so h{x) G dY . Therefore, 
dist(y,9Y) = dist{hj (x) , dY) ^ \hj{x) - h{x)\ 

^ \\hj - ^||^°°(K) ^ sup \\hj - /i||jfoo(K) = ej. 

Hence y eYj, as claimed. 

Now it follows that 

/ Jh{x)dx ^ lim \Yj\ 
Jk i^'^ 

because we have a decreasing sequence Yi D Y2 D . . . of measurable subsets 
of Y whose intersection is empty. This yields 



Jh{x) dx = 0. 



Finaly, B can be expressed as a union of increasing sequence IK' C IK" C 
K'" C . . . of compact subsets of B and a set E C i3 of measure zero. We 
see that 



Jh= Jh+ A = + = 0. 

IB JE Jk'uk"u... 
Since Jh ^ almost every in B, wc conclude Jh = in B. □ 

6.0.1. The 6 -function. Let X and Y be bounded doubly connected domains. 
We recall the notation X/ and Xq for the bounded (inner) and unbounded 
(outer) components of C \ X. Thus dXj and OXq are exactly the boundary 
components of X, so C = X/ U X U Xq- The same notation applies to Y. 
Let h e J^i^^^ (X, Y). We define a function 5:C°-^[0,1] by the rule 

'0, ifxeXi 

(6.2) S{x) = { dist[/i(a;) , aVj] Vdist[fe(x) , OYq] ' X E X 

1 , if a; G Xo . 

Lemma 6.2. The S-function is continuous and monotone. 

Proof. The mapping h G J^^;(X,Y) is a weak limit of homeomorphisms 
hj : X ^ Y in the Sobolev space W^''^{X , Y). Wc appeal to [28, Theorem 
1.1] which provides us with a uniform bound of the distance function of any 
homeomorphism /: X Y in the Sobolev space W^''^{X , Y) ; namely, 

dist [/(x) , dY] ^ T]^^{x) <%[/] for all xeC 

where 77^^ : C ^ [0 , 00) is (independent of / ) a continuous function van- 
ishing outside X . Thus we have well defined continuous monotone functions 

'0, ifxeXi 

X.(^\-} dist[hj{x),dYi] ifT-f=X 

"jK-^J — S dist[hj{x),dYj]+dist[hj{x),dYo] ' t ^ 

1 , if x G Xq ■ 



DIRICHLET ENERGY AND HOPE DIFFERENTIALS 



23 



converging uniformly to 6{x). Therefore, the 6 -function is also continuous 
and monotone. □ 

This implies that the set h~^{Y) = {x e C: < S{x) < 1} is connected 
and its complement consists of two components 

C \ h-\Y) = {x: S{x) = 0} U {x: S{x) = 1} . 

These arc disjoint continua containing d%j and SXq, respectively. We have 
established the following fact. 

Lemma 6.3. Let X and Y be bounded doubly connected domains and h G 
Mll^iX,Y). Then h-\Y) is a doubly connected domain separating the 
boundary components o/X. 

Proposition 6.4. Every energy-minimal map h G J^2(X, Y) is a harmonic 
diffeomorphism of h'^^M) C X onto Y. 

Proof. This is immediate from Proposition 4.1 on harmonic replacements. 

Indeed, we may cover X by inner cells. In each inner cell U C X the 
map /i is a diffeomorphisms, exactly equal to its harmonic replacement, 
since otherwise we would have a map in J^2(X, Y) with smaller energy. It 
remains to observe that a local diffeomorphism which is a c - uniform limit 
of homeomorphisms must be a global diffeomorphisms. □ 

The following existence and uniqueness result is from [20, Theorems 2.3 
and 2.4] 

Theorem 6.5. Let X and Y be doubly connected domains so that Mod Y ^ 
ModX. Then there exists a '^'^-diffeomorphism H:X°-^Y in W^^'^{X,Y) 
that minimizes the Dirichlet energy among all mappings in S(X, Y). More- 
over, such a minimizer is unique up to a conformal change of variables in 
X . //, however, Mod Y is sufficiently small relative to Mod X , then no 
-diffeomorphism of X onto Y is energy-minimal. 

7. Cracks in doubly connected domains 

Let X and Y be bounded doubly connected domains of finite conformal 
modulus. We shall establish basic geometric features of the solutions to the 
Hopf-Laplace equation 

(7.1) hJT^ = (p for /i G >r^'2(X, Y) 

where (p G =Sf^(X) is analytic. Without any additional assumptions of topo- 
logical nature, this equation admits rather bizarre solutions, see §3.6. Our 
main interest, however, is in the energy-minimal solutions. Prom this point 
of view it is natural to restrict the equation (7.1) to the class =^i^^(X, Y) of 
weak limits of homeomorphisms hj : X Y in the Sobolev space y^^'^(X, Y) 
We also assume that the Hopf differential hzh^dz <S> dz is real and nonva- 
nishing along dX. The simplest way to rigorously formulate this property 
is through a conformal transformation / : X A of X onto an annulus 
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A = {z: r < \z\ < R}. Upon such transformation the Hopf differential 
remains real along the boundary of A , by definition. For X = A this reads 

as 

(7.2) hjiz = ^, cGM, c^O. 

Z'^ 

Thus the rays 71? = {pe*^ : r < p < R} , O<0<27r, and the concentric 
circles Cp = {pe'^ : ^ 9 < 27r}, r < p< R, are the orthogonal families of 
the trajectories of h^hzdz ® dz. Their preimages under /: X A, also 
known as isothermal curves, are the trajectories in X. Note a dissonance 
that the differential 

z^ 

is positive along 5X if c < and negative along OX if c > . 

Theorem 7.1. Suppose that the differential hzhz dz^dz for h G J^^^ {X,Y) 
is negative along the boundary o/X. Then h: X Y is a harmonic dif- 
feomorphism. 

Proof. Applying a conformal change of variables we may assume that X is 
an annulus A = {z: r < \z\ < R}. Upon such a transformation the Hopf- 
Laplace equation reads as 

c 

hzhz = for some c > 0. 

First we show that h takes A onto Y. For this, recall the inclusions Y C 
h{A) C Y, see (6.1). Suppose to the contrary that there is Xo G A such that 
= yo := h{xo) G dY. We invoke the J-function given by (6.2). Accordingly, 
5{xo) is either equal to or 1; say 5{xo) = 0. We look at the zero set 5~^(0) = 
{x G C: 5{x) = 0}. Clearly, it contains A/ and the point Xo, r < \xo\ < R. 
By Lemma 6.2 the function 5 is continuous and monotone, so (5~^(0) is a 
continuum in C. The circles Cr = {z: \z\ = r}, r < < |xo| , separate 
Xo from A/, thus intersect 5^^(0). These circles arc vertical trajectories 
of ^p{z) = cz^'^ , c > 0. We now appeal to Proposition 3.11 which tells us 
that the union Ur<r<|a;o| measure, a clear contradiction. Thus 

h{K) = Y . Finally, we apply Theorem 1.12 in [14] which asserts: 

Theorem 7.2 (Partial Harmonicity) . Every deformation h G S(X, Y) be- 
tween bounded multiply connected domains, which satisfies the Hopf-Laplace 
eqyation is a harmonic diffeomorphism of h^^{Y) C X onto Y . 

We then conclude that h : A Y is a harmonic diffeomorphism. □ 

The situation is dramatically different if hzhz dz(Sidz is positive along 9X; 
cracks along vertical trajectories emanating from the boundary components 
and ending inside X are likely to occur. A complete description of such 
cracks is given in Theorem 7.3 below where we assume, without loss of 
generality, that X is an annulus. The rays and concentric circles of the 
annulus are vertical and horizontal trajectories. 
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Theorem 7.3. Suppose that h G J^^^{A,Y) satisfies the Hopf-Laplace 
equation 

c 

hzhz = for some c < 

in an annulus K = {z: r <\z\ < K]. Then every ray TZ^ = {pe*^: r < p < 
R}, ^ < 2tt , consists of three disjoint subintervals (some can be empty) 







7^^ = 7^? 






= {pe'' 


r < ps^re] 


inner ray (possibly empty) 




= {pe'' 


rg < p < Rg} 


middle ray 


K 


= {pe'' 


Re^P<R} 


outer ray (possibly empty) 



where r ^ rg < Rg ^ R. The union of middle intervals, denoted by 

A= U n^A 

Os;6»<27r 

is a doubly connected domain separating the boundary circles of A. More- 
over, 

(i) h: A Y is a harmonic homeomorphism 
(a) h: IZ^j — )• dYi is constant (if not empty) 
(Hi) h: TZq —7- dYo is constant (if not empty) 




Figure 1. Cracks in a doubly connected domain and their 
conformal images (rays) in an annulus. 



In what follows we denote the doubly connected domain /i^^(Y) by A. 
Furthermore, Aj stands for the bounded component of C \ ^ and Ao for 
the unbounded component of C \ ^. We write B := A \ Br := Aj Ci B 
and Br := Aq H B. In view of (6.1) B = h'^{dY). By Lemma 6.3 the sets 
Br and Br are disjoint and relatively closed in A. 



26 



T. IWANIEC AND J. ONNINEN 



Remark 7.4. Knowing h: A Y will allow us to determine uniquely the 
mapping ^: A — >■ Y. 

The proof of Theorem 7.3 will employ the following lemmas. 

— 1 2 

Lemma 7.5. If h: A ^ Y belongs to ^xxm (''^j^) '^'^'^ V° ^ then the set 
h~^{yo) U A/ is connected. Similarly, if yo G Yq, then the set h~^{yo) U Aq 
is connected. 

Proof. The problem translates into a question about monotone mappings 
of §^ onto itself. We choose and fix homeomorphisms tp: A and 
i/j: Y where is the 2-sphere with two points removed, say the 

north and south poles. We assume that ip takes dAj and dAp into north 
and south poles, respectively, in the sense of cluster limits. Similarly, we 
assume that ^ takes dYj and OYq into north and south poles, respectively. 
Every homeomorphism / G J^(A,Y) induces a unique homeomorphism 
F : §2 §2 g^(,j^ ^j^g^^ F o = tp o f. Note that F takes the north pole 
into the north pole and the south pole into the south pole. Moreover, if a 
sequence of homeomorphisms hj G J^j^^(A, Y) converges weakly to h, then 
upon lifting to we obtain homeomorphisms Hj-. converging 
uniformly to a surjective mapping H : S^. This is immediate from 

the uniform bounds 



dist[hj{x),dYi] ^ ri^{x)^^x[hj] 
dist [hj{x),dYo] ^ Voi^)\f^\hj\ 



and 



\hj{xi) - hj{x2)\ ^ 



2 ^ C'k 



g _j_ diam K 



log (t-^i^^^ 

for xi,X2 G IK, where K is any compact subset of A. Here the functions 
rjj , tJq G "^(M^) are nonnegative, r]j{x) = on A/ and r]Q{x) = on 
Ao , see Lemma 2.9. The limit map H : S keeps the north and south 

poles fixed and is monotone, see the Kuratowski-Lacher Theorem 2.8. In 
particular, if C is a continuum in then so is H^^{C). Now Lemma 7.5 
follows if we return to the domains A and Y via the inverse homeomorphisms 
if-^ : ^ A and ip-^ : ^ Y. □ 

Lemma 7.6. Suppose that h G J^^^{A,Y) and 

hzhz = for some c < 0. 

If C C A is a nonempty connected set such that h\^ = {yo}- Then C C 
IZ^ = {pe'^ : r < p < R} , for some 9 <2tt. 

Proof. Consider the projection p: A — ^ p[z) = |||. Thus p{C), being a 

continuous image of a connected set, is a connected subset of S^. If p{C) is a 
single point say e'^, for some 6, then C dTlP . Suppose to the contrary that 



DIRICHLET ENERGY AND HOPE DIFFERENTIALS 



27 



there is an arc {e*^: a < 6 < P} C p{C). Thus every ray TZ^, a < 6 < P, 
intersects C. On the other hand, by Proposition 3.11, the union of rays 
which intersect h~^{yo) has zero measure, a contradiction. □ 

Proof of Theorem 7.3. Consider a ray TZ^ C A. There is a point z = pe^^ G 
A, for some p G (r, R), since otherwise the set = AiUAq^TZ^ would be 
connected, contradicting double connectivity of A, see Lemma 6.3. Denote 

re = mi{p G (r, R) : pe'^ G A} 
Re = sup{/9 G (r, i?) : pe'^ G .4} 

so we have r ^ rg < Rg ^ R. Associated with these numbers is the partition 
of TZ^ into: left-open right-closed interval TZ^j = {pe*^ : r < p ^ vg} (possibly 
empty), the open non-empty interval TZj^ = {pe*^ : rg < p < Rg}, and left- 
closed right-open interval TZq = {pe*^ : Rp ^ p < R} (possibly empty). 

Claim. Every Zo = PoC*^ G TZ^^ ; that is, rg < po < Rg , lies in A. 
Suppose to the contrary that h{z) G dY, say h{z) = yo G dYo (the case 
yo G dYi is similar). Consider the set 

C = h-\yo)cA 

and its all connected components {C,y},y^Mi so C = \J Ci,. Let Co 

denote the component that contains Zo- According to Lemma 7.5 the set Co 
intersects dAp- On the other hand, by Lemma 7.6, Co lies in one and only 
one of the rays of A, so in TZ^. Thus Co , being a connected set, takes the 
form 

Co = {pe'^ : R'e^p<R} for some r < R'g < R. 
Since poe'^ G Co, we see that 
(7.3) R'e^poKRg. 

On the other hand, since /i(pe*^) = yo for every Rg ^ p < R and h{pe^^) ^ yo 
for every rg < p < Rg, it follows that Rg ^ Rg, contradicting (7.3). The 
claim is verified. 

Now, by the definition of Rg, no point pe*^ with Rg < p < R lies in A. 
Thus the interval TZq lies inA\^ = ^ = ;S,.U Br. Actually TZq has to 

lie in the outer component of C \ .A, the point i?e*^ belongs to TZq . Thus 
C A n and h{pe'^) = yo = h{Ree^^) for every Rg ^ p < R. In 
the same way we see that TZ^j C AD Br and h{pe^^) = h{ree^^) for every 
r < p ^rg. □ 

8. An integral identity 

Lemma 8.1. Let X, Y and G be bounded domains in C. Suppose that 
h: G Y andH: X Y are orientation preserving -diffeomorphisms 
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of finite energy. Define f = o h : G X. Then we have 

'\f.-iiz)M' 



-S'G[h]=4: 
(8.1) ''-'^ 



l/.P-l/.-P 
{\hz \ - \h,\f- \fr^' 



\hzhz\ dz 



+ 4 // ^'-^"....'-^"/..o"" dz 

J JG 



where 

\hzhs\hzhz\~^ if hzhz 7^0 

1 otherwise. 



7 = 



The integrals in (8.1) converge . 

Proof. It is worth noting that / : G X need not have finite energy. The 
convergence of the integrals, not obvious at the first glance, is a consequence 
of the finite energy condition imposed on the mappings h and H . 
We begin with the chain rule applied to H = h o f~^ : X Y , 



dHjw) ^ j^j.dr^ ^ h-Jz)^ 
dw ^ dw ^ dw 



dHjw) ^ ^j.dr^ + h-Az)^ 
dw ^ dw ^ dw 

where w = f{z) . We express the complex partial derivatives of f~^: X ^ X 
at w in terms fz{z) and fz{z) at z = f^^{w) , 

df-' _ JJV) df-^ _ fz{z) 



dw J{zJ) dw J{z,f) 

Note that the Jacobian determinant J{z,f) = — l/^p is strictly posi- 
tive. These expressions yield 

dH hzfz - hzfz , dH hzfz - hzfz 
and 



dw dw |/.|2-|M2 

Next we compute the energy of H over the set /(G) = X by substitution 

w = x(.z), 

^f(^G)[H] = 2 // {\H^f + |i?^|2) dw 
JJfiG) 

2 ff \hzJz - hzfz? + \h-zfz - hzfzl^ 

J JG 



On the other hand, the energy of h over the set G equals 



SG[h] =2 ff {\hz? + \hzf) dz 
J JG 



DIRICHLET ENERGY AND HOPE DIFFERENTIALS 



29 



The desired formula follows by subtracting these two integrals, 



S.[H] - SG[h] = Ajj^ ^^-^^p d. 

/X \fz?-\h? 



^ I I 2\hzhz\ ■ l/zp - 2 Re [hzhzfzXz] 



(8-2) +4 [[ „ ,J"^'/„ dz 

JJG 



4 



G \fz?-\h? 

\f.-l{z)h? 



1 



l/.p-l/.-p 

G \fz?-\h? 



hzhzl dz 



+4 // d^ 



□ 

9. Some Free Lagrangians, Normal and Tangential Distortions 

In the 1980-ies a novel approach towards minimization of polyconvex en- 
ergy functionals for mappings between domains in was developed and 
published by Ball [6]. The underlying idea was to view the integrand as 
convex function of null Lagrangians. The term null Lagrangian pertains 
to a nonlinear differential expression whose integral over any open region 
depends only on the boundary values of the mapping, sec [10, 16, 18]. In 
this paper the mappings of our interest are homeomorphisms h: X Y 
that are given only on a part of the boundary (possibly empty part, in 
which case we are dealing with a traction free problem). There exist some 
nonlinear differential forms, defined on a class of homeomorphisms, whose 
integral means remain independent of their boundary values. These are 
rather special null Lagrangians, called free Lagrangians [27]. Precisely, a 
free Lagrangian for a pair of domains X, Y C M" is a nonlinear differential 
n -form L{x,h, Dh) dx whose integral mean over X depends only on the 
homotopy class of a homeomorphism h: X Y. Here are a few of them 
in the planar domains, where we make use of polar coordinates p and 9 

(9.1) z = pe'^ , ^ p < GO and ^ 6* < 27r . 

The normal (radial) and tangential (angular) derivatives of a Sobolev map- 
ping / are defined by 

dfipe'^) 



(9.2) fNiz) 

up 

and 

(9.3) f^^'^ -p^e- ^ ^=l"l- 

The Jacobian determinant of / is 
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The following three free-Lagrangians are concerned with a pair of the annuli 
A = {x G C: r < < i?} and A* = {y G C: r* < |y| < i?*}. 

• Fullback of a 2-form in Y via an orientation preserving homeomor- 
phism h € J^{X,Y) n ^^^'^(X, Y) is a free Lagrangian 

(9.4) / / N{\h\) J{x, h)dx = [[ N{\y\) dy . 

• Normal differentiation gives rise to a free Lagrangian for h G ^(X, A*)n 
>ri'i(X,A*) defined by 



(9.5) 



\h\ 



N 



\X\ 



dx 



27r 
27r 



f.(|.,)fd. 



A{t) dr 



(9.6) 



• A dual free Lagrangian for h G (A, Y) fl >^-'^'^(A,Y) arises from 
tangential differentiation 

■R 

r \J\x\=t 
■R 

27r / B{t)dt 



dAxgh 

de 



de dt 



Lemma 9.1. Let H be a bounded doubly connected domain that separates 
the origin from oo , and let A* = 74(r*,i?*) be a circular annulus. If 
he je(X,A*)nW^'^{X,A*) then 

—r^ ^ 27rlog(i?*/r*). 

Proof. Choosing A{t) = l/r in (9.5) we have 

^R, dr 



IIx 



27rlog(i?*/r*) = 27r 



and 



\hN\ dz 
\h\ \z\ 

|/lAr| dz 



ff\hN\^y < ffiM^ ff = ff\M^ ff 

JJx \h\ \z\J ^Jh Jh kP77x|/iP JJx Jh \z\'JJa^ 



\y\' 



by (9.4). 



□ 



Lemma 9.2. Let A = A{r, R) be a circular annulus and Y a bounded doubly 
connected domain of finite conformal modulus. Ifh G ,;?^(A, Y)n>^^'^(A, Y), 
then 



IIa 



[MP ^ \og\R/r) 
A Jh \z\^ ^ ^ ModY ■ 
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onto. 



Proof. There exists a conformal transformation F: Y A* , of Y onto 
an annulus A* = {z: < r* < |z| < i?*}. We define g = F o h: Y 
Since F is conformal 



rr \hT\' dz _ rr 
JJk Jh JJh 



J9 I^P 



Applying (9.6) with B{t) = 1/t we have 

IIa 



dt 

27rlogi?/r = 27r ' 



Now, it follows by Holders inequality that 



ViA 1^1 \z\J ^ JJa Jg iA I5P JJa Jg iA* 



□ 



Lemma 9.3. Let A = yl.(r, R) be a circular annulus, < r < R < oo, and 
Y a bounded doubly connected domain. Suppose h G >^^^^(A,Y) satisfies 
the Hopf-Laplace equation 

Q 

(9.7) hzhz = in A , and J/j ^ almost everywhere 

where c G M is a constant. Then we have point-wise inequalities 

(9.8) 



Proof. The complex Hopf-Laplace equation (9.7) reduces to the system of 
two real equations, 

(9.9) I^ATp - Ihrf = ^ ; 

(9.10) Re(h^hT) = 0. 

Recall that Jh = Im/ijv /it ^ which in view of (9.10) reads as 

(9.11) Jh = \hN\\hT\ 

Combining this and (9.9) the inequalities (9.8) follow. □ 
Geometric function theory is concerned with the distortion function 

" = — = J, • 

For the purpose of this paper we will decompose it as = + , 
where (using polar coordinates) the normal and tangential distortions of / 
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are defined by the rules 



(9.12) K' ■ '^"^ 



(9.13) k{ ■= ^— = ^ 



^/■-|2 lf^|2 



By convention, these two quotients are understood as whenever the nu- 
merator vanishes. Naturally, they assume the value +cx) if the Jacobian 
vanishes but the numerator does not. For a mapping / G the quanti- 

ties /tv , /t J and J/ are finite a.e. and, therefore, K^^ and Kj, are well 
defined measurable functions in the domain of definition of / . 

9.1. The distortion of the difference of two solutions. Suppose /i, G 
>^^'^(X,C), J/i ^ and Jh ^ , have the same Hopf-product, 

hzhz = HzHz = ip{z) ^ almost everywhere (not necessarily analytic). 

Consider the difference 

F{z) = H{z) = h{z) G /r^'2(X,C). 

We have 

hzFz = hz{Hz — hz) = hgHz — ^ = hzHz — HzHz = Hz{hz — Hz) = —FzHz 
where we notice that 

< \hz\\hz\ = \ip\ 
\H-z\^ ^ \Hz\\H-z\ = \^\ 

Hence ^ li^plF^p so Jj? ^ almost everywhere. Next we introduce 

the Beltrami distortion coefficients 

^h{z) = ^ 1 and kniz) = ^ 1. 
\nz\ l-Wzl 

We find that 



\Fz\ = kF{z)\Fz\ where kpiz) = \/kh{z)kH{z) ^ 1. 
Indeed, we have 

\hz\'^\Fz\'^ = \Hzf\Fzf , where l^^-p = and kaWzf = \^\. 

Hence 

khkH\Fzf\v\ = kH\hz\^\Fz\^ = V^l' = V \Fzf 

and therefore 

\Fz\ = \fkhkH\Fz\ = kF{z)\Fz\. 

Note that 

< 1 whenever Jh + Jh ^ 
= 1 whenever + Jh = 0. 



kpiz) 
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In particular, F has finite distortion whenever Jh or Jh{z) ^ 

\DF\^ = 2{\Fz\^ + \F,\^) = 2\±^Jp . 
^ ' l-khkn 

10. HOPF SOLUTIONS ARE ENERGY-MINIMAL, PROOF OF THEOREM 1.4 

Throughout this section X and Y are arbitrary bounded doubly connected 
domains, X being nondegenerate. We only need to prove the "if " part of 
Theorem 1.4. The "only if " part follows from Corollary 3.6 , because the 
energy-minimal mappings are critical points for all variations along 3X . 
In the case of negative Hopf-differential we actually obtain slightly better 
result; it holds in a larger class of the Hopf solutions. Precisely, 

Proposition 10.1. Suppose the Hopf-differential hzh^dz dz defined for 
h G Y) is holomorphic and real negative along dX. Then 

(10.1) <^x[h] = inf{ ^xb] : 9 G (X, Y) } . 

Furthermore, h is a unique ( up to the conformal change of variables inX) 
minimizer in ^j^j^ (X, Y) . 

Obviously, if h happen to belong to !^2(X,Y) C Ji^l^{X,Y), then 
the equation (10.1) yields (1.3), as desired. However, as yet, the equal- 

1 2 

ity J^ij^(X,Y) = J^2(X,Y) has been established only under Lipschitz 
regularity of Y, [29]. 

Proof. By virtue of Theorem 7.1 h: X Y is a harmonic diffeomorphism. 

A conformal transformation of X onto an annulus A = {z: r < \z\ < R} 
takes ^jj^ (X, Y) onto (A,Y). Thus it involves no loss of generality in 
assuming that X = A, so we have the equality 

hzhz = ^ c G M \ {0} . 

The assumption that Lp[z) dz^dz is negative along dA yields c > 0. In view 
of Lemmas 9.3 and 9.2 we have 

„ log2 T ff I^tP dz ff dz ^ ^ R 

This shows that Mod A ^ ModY. Now we appeal to Theorem 6.5 which 
asserts that there exists a C°°-diffeomorphism H : A '^^^ Y such that 

^x[H] = inf{ ^x[g] : <7 G (X, Y) } . 

Furthermore, H is unique (up to the conformal change of variables in X) 
minimizer in J^^^^ (X,Y). It remains to show that SxlH] = S^ih]. For this 
purpose we consider a diffeomorphism 

f:=H-^oh: A°-^A. 
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Applying Lemma 8.1, we have 



_ 2 f_\2 

<rA[i^] -<^A[/i] =4c ~ 

(10.2) 



\fz 2 1 ^ 



l/.P-l/.-P 



dz 

lTi2 



2 I -p |2 



{\h,\ - \h,\Y\f. 



+ 4 dz 



l/.P-l/.-P 
Since 

|/z - f /gP _ / 

l/.p-l/.-p"^ 

by Lemma 9.2 we conclude with the inequality 

l/.-f/.-p dz ^ , , , /•/■ dz 



UP' 



Furthermore, / is an orientation-preserving diffeomorphism. Therefore, 
(10.3) Sk[H] - > 0, 

and so Sk[H] = S^M- □ 

It is the case of positive Hopf-differential that, for the proof of Theo- 
rem 1.4, we need to take the solutions in the class of strong limits of home- 
omorphisms. 

Proposition 10.2. Suppose that a Hopf-differential hzhzdz ® dz , defined 
for h E 2(X, Y) , is holomorphic and real positive along dX . Then h is 
an energy-minimal. 

Proof. A conformal transformation of X onto an annulus A = {z: r < \z\ < 
R} induces an isometry of J^2{'X, Y) onto Jif2{A, Y). Thus we may assume 
that X = A, so as to apply the equality 

hzhz = ^ c G M \ {0} . 

The assumption that (^(z) dz ® dz is positive along dA simply means that 
c < 0. _ 
We write G := ^-^(Y). In view of Theorem 7.2 the mapping h: A Y 
is a harmonic diffeomorphism from G C A onto Y. Let H : A Y be an 
orientation preserving ^°°-diffeomorphism. We denote 



f = H-'^oh:G°-^A. 
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l/z + f /z| 
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SK[H\-SG[h]=4.\c\ 



Ll/.P-l/.-P 



1 



dz 



(10.4) 



+ 4 
= 4 Id 



G 



K-^ — 1 



dz 



+ 4 



IIg 



Jf 



dz. 



Before estimating the right hand side we will show that 

dz 

M2- 



J Ja\g 



(10.5) 

In view of Lemma 6.1 J/j = in A\ G. Since c < 0, by Lemma 9.3, ]/7.jv| = 
in A \ G. Therefore, = |/iTp in A \ G. On the other hand, in view 

of (9.9), we have |^tP = -4c|z|-2 in A\G. Therefore 

// \Dh\' = -Acf A. 
JJa\g Ja\g \z\ 

Combining (10.4) with (10.5) we arrive at the identity 



^a[H] - 



(10.6) 



= 4|c 
+ 4 



Ki{z)dz^ 
I - \h,\?\h 



Jf 



dz 

A I^^P 



dz. 



According to Lemma 9.1 



// i^^^27rlog 
J JG 



Therefore, if ii": A — 
(10.7) <^a[H] - ^aW ^ 4 



Y is a ^°°-diffeomorphism, we can write 

iM-\h,\?\h ^' 
G Jf 



dz>0. 



The last inequality follows from the fact that / preserves the orientation. 
Hence S'aIH] ^ <^a[^] for an arbitrary H G J^(A,Y), meaning that h is an 
energy-minimal map. 

11. Uniqueness, Proof of Theorem 1.3 

Let H and h be energy-minimal mappings in (X, Y) . According to 
Theorem 1.4 the Hopf differentials of h and H are analytic and real along 
dX. 

Case 1 If the mapping h has negative Hopf-differential along , then 
the uniqueness ( h differs from if by a conformal change of variables in X ) 
was already established in Proposition 10.1. 
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Case 2 Suppose the mapping h has positive Hopf-differential along dX . 
We may assume that X is an annulus. Therefore, 

Q 

hzhz = for some c < 

c 

HzHz = -|- for some G M. 

The reader may wish to notice that at this point we do not know (though 
it win be shown to be true) that the real valued constant coincides with 
c. We also cannot assume that < 0. Let us denote by G := h~^{Y) 
and := H~^{Y). These are doubly connected domains separating the 
boundary circles, Lemma 6.3. In view of Theorem 7.2 the mapping H 
is a harmonic diffeomorphism from onto Y. Recall that G = h~^{Y). 
Thus / := o h: G Gh is an orientation preserving diffeomorphism 
which also preserves the order of the boundary components. We denote 
the inverse of f hy g = = o H: Gh G. Choose and fix 
a disk D C B C Gh ■ The aim is to prove that g is analytic in D and, 
therefore, analytic in Gh ■ Since H G 2(A, Y) = Diff^ (A, Y) there exists 
a sequence of diffeomorphisms ■ A '^^^ Y, k = 1,2, . . . , converging to H 
in #'^'^(A) and c-uniformly. In analogy to / and g we define 

f'' ■= i^-i oh ■.Gh°-^A and g'' := H'^ oh :A°-^Gh. 

Since iJ^: D ^ Y converge uniformly to H : B H{B) , where H{B) is 
a compact subset of Y , there is a neighborhood U of iJ(D) , compactly 
contained in Y, such that i?jfc(B) C U for all sufficiently large k, say 
k^ ko- Since U is compact in Y the set F := /i~^(U) is compact in G. 
Then we note that 

/(D) = h-^{Hk{^)) C /i~^(U) = F for /c ^ K. 

Furthermore g^ converge uniformly to g = h'^ o H : ^ h'^ {H{Ii)) C F . 
In view of (10.7) we have 

JJg '^f' 

J Jw Jf' 

r r I ffe|2 

^ 4g2 / / \Jij_ g > 0. 

J Jv "J 

Here we used the inequality 

{\hz\ - \hz\f ^ > for z G F 

because h is a, difi'eomorphism on G . In the later integral we will make a 
substitution z = g^{w) , where w G /^(F) , to obtain 

^A[Hk] - ^A[h] >^e^ If \9U^)\'' dw > 1 1 \g%{w)\'' dw 
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Letting /c — )• oo we find that ^ in (D) . Since ^ g uniformly 
on D we see that = on D , where we recall that B was an arbitrary 
open disk compactly contained in Gh- Consequently, g: Gh G is a 
conformal map, so is / : G Gh ■ This conformal map preserves the 
order of the boundary components, so 



z)dz 



27riJ^ f{z) 

where 7 can be any smooth Jordan curve that separates the boundary com- 
ponents of G, and is oriented counterclockwise. Next wc differentiate the 
equation = h{z) with respect to z G G . Chain rule yields, 

dH ^ , _ , dH- 



^^-f{z) = K{z) and -jf'{z) = h,{z). 

These formulas together with the following two Hopf-Laplace equations 
dHdH c„ ^ ^ 7- c 

yield a differential equation for /, 

(zf\^ c 

(11.1) I — — = — ; this is a nonzero constant. 

V / y c„ 

Therefore, we have a continuous branch of the square root 
— ^ = a a constant (real or imaginary). 
On the other hand, 

_ 1 r adz _ I rf'{z)dz_^ 
2m z 2TTi f{z) 

Therefore, equation (11.1) reduces to = in G. Thus f = Xz for some 

constant A / 0. Substituting into (11.1) we find that c = c^. We may, and 
do, assume that A is a real positive number, simply by rotating the z-variable 
H = H{z) if necessary. We may also assume that < A ^ 1, because 
H{z) = h{z/X) for z G G^ and the mappings H and h are interchangeable. 
Summarizing 

G^ = AG C A, < A ^ 1, H{Xz) = h{z) for z G G and c = c^. 
Now come the key ingredients of the proof. Consider the difference 

F{z) = H{z) - h{z). 
Claim. For every ^ 6 < lir we have 

(11.2) lim F{pe'^) = = lim F{pe'^) 

p\,r p /■r 
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Proof of Claim. We recall that 

H{Xz) = h{z) for z G G 

and = AG C A for some < A ^ 1. Using the same notation as in 
Theorem 7.3 we recall the ray segments 'RP and the corresponding partitions 

r ^ rg < Ro ^ R for ^ 
r^Xrg<XRe^R for iJ. 

We shall only demonstrate the proof of the limit in (11.2) at the inner 
boundary, the proof for the outer boundary is similar. 

The case A = 1 is obvious, because G = Gh and H{z) = h{z) for z e G. 
According to Remark 7.4 the values of H and h outside G arc uniquely 
determined by those on G, so H = h on A. Now assume that < A < 1. 

Case 1. Let r < Ar^* , so the intervals {r,rg] and (r, Xrg] are not empty. 
We have, for every r < p ^ Xro < ro, the following equalities, 

H{pe'^) = H{Xrg e'^) = h{ree'^) = h{Xrg e'^). 

Therefore, 

F{p e'^) = for all r < p ^ A rg. 
Case 2. Let r = Xrg < rg. For every r < S ^ ro we have 

h{d e^^) = h{re e*^) = lim h{S e^^), because re e'^ eGn A. 

S\re 

Since 5 e'^ G G we it follows that 

h{ree'^) = lim h{XSe'^) = lim H{pe'^). 

S\re p\r 



Thus 



limF(pe^^) = 0. 



□ 



The function F = H — h & #'^'^(A), admits an extension to the entire 
complex plane as a member of >^^'^(C). In particular, F restricted to almost 
every ray 

TZ^ = {pe'^: r < p < R}, 6 < 2tt 

has an extension as an absolutely continuous function on TZl = {p e*^ : r* ^ 
p ^ R*} , where 0< r*<r<i?< J?* <oo. Since 

lim F{p e*^) = = lim F{p e^^) 

by (11.2) the zero extension of -F: A C to the entire complex plane, 
denoted by Fo : C ^ C , is absolutely continuous on almost every ray TZ^ = 
{pe'^ : ^ p ^ R} . Obviously, Fo is absolutely continuous on almost every 
circle Cp = {p e*^ : ^ 61 < 27r} . In other words Fo G W^''^{C) and vanishes 
outside A . This yields 



// Jf= [[ Jf.=0. 

JJA J JC 
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Since Jp ^ almost everywhere in A we conclude that Jp = almost 
everywhere in A . Next for z eGU we have the distortion inequality 

Thus DF = almost everywhere on G U which implies that F = 
constant on every component of G U G^ . We claim that G = G^ . Suppose 
to the contrary that G^ \ G 7^ , so there is a disk D C G^ \ G , and 
F = constant on D . Let 7^ C D be a ray segment. Since D C A \ G 
the function h is constant on 7^ . In particular, H = F + h is constant 
on I C G^ , contrary to the injectivity of H when restricted to G^ . In 
summary, G = G^ , so A = 1 and H{z) = h{z) on G . This yields H = /;, 
in A , by Remark 7.4. □ 

12. How TO MAKE A SQUARE HOLE WASHER FROM AN ANNULUS 

To illustrate the ideas developed above we shall take on stage the energy 
minimal deformations of an annulus 

^ = {xi + 1X2-. l< x\ + xl < R^} 

onto a square washer 

Y = {yi + iy2: 1< + |y2|, vl + vl < P^} , P>1 

The method of exploring symmetries in this example actually applies to 
hexagonal nuts or many more symmetric doubly connected polygonal figures. 
We shall look for the energy-minimal deformations among strong W^'^- 
limits of orientation preserving homeomorphisms h : X Y . Recall the 
notation ^ (X , Y) for the class of such mappings. By virtue of Theorems 
1.1 and 1.3 , we see that 

Theorem 12.1. There is a map H G J^(X,Y) (unique up to the rotation 
of X) with least energy. 

The geometric features of the energy-minimal mapping h depend on the 
outer radius of the annulus X , so we vary it while keeping the target Y 
unchanged. 

Theorem 12.2. There are two critical outer radii 1 < Ro ^ R° such that: 

1. For 1 < R ^ Ro , the energy-minimal map h is a harmonic 
dijfeomorphism of X onto Y ; no cracks occur in X . 

2. For R^ R° , h fails to be injective; cracks are unavoidable. 
For all R > R° , the energy-minimal deformation has the following 
properties: 

3. There are four straight line segments (possibly points in dXj ) 
I, il, I, —il , where /=[1, a]cM ,1 ^ a < R , such that 

• h{T) = {1}, h{il} = {i}, h{-I) = {-1} and h{-il) = {-?;} 

• The set X' = X \ [ / U i/ U — / U —il\ is a doubly connected 
domain and h : X' "^^^ Y is a harmonic diffeomorphism. 
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• The map h is invariant under the rotation by the right angle; 
h{iz) = ih{z) for all z € X . 

!/2 




Figure 2. How to make a square hole washer from an annulus. 

An interpretation of this result is that in order to make a square hole 
washer from a thick annulus one must make cuts in the annulus. The 
energy-minimal deformation takes those cracks into the corners of the inner 
boundary of Y , exactly where it fails to be convex. 

Theorem 12.3. In spite of the presence of cracks in case R ^ R° , the 
energy-minimal map is Lipschitz continuous up to the closure X . 

12.1. Proof of Theorem 12.2. 

Proof. Conditions 1 and 2 are immediate from Theorem 6.5. In the proof 
of Condition 3 we shall explore the symmetries in both the domain X and 
the target Y. We shall first construct an energy-minimal harmonic diffeo- 
morphism of the first quadrant sector 

S = {pe'^ -.KpKR, 0<e < tt/2} 

onto the first quadrant of the square hole washer, 

T = { yi + i y2 : 1 < yi + 2/2 , yl + yl< , yi > 0, 2/2 > } . 

We view S and T as quadrilaterals with corners. 

1, i?, ii?, i G and 1, i G OT, respectively. 

Consider the class Q2(S,T) of all quadrilateral mappings from S onto 
T . These are homeomorphisms / : S T which extend continuously to 
/ : S T and take the corners of S into the corners of T in the following 
order 

/(1) = 1, f{R)=p. fiiR) = ip, f{i) = i- 
Lemma 12.4. There exists h G Q2(S,T) of smallest energy. 
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Remark 12.5. The interested reader may wish to notice that this lemma and 
its proof given below actually hold for every pair of quadrilaterals in which 
the target is convex. 

Proof of Lemma 12.4. Let {/j}j=i,2,... denote the minimizing sequence of 
quadrilateral mappings fj : S T , so 

^ inf{ 4[/] : / e Q2(S,T) }^m. 

We assume that {fj} converges weakly in to some mapping / G 

y^^'2(S) and observe that fj : S T also converge uniformly to /: S ^ 
T , see Remark 2.10. 

The boundary map / : dS dT , being a uniform limit of homeomor- 
phisms fj : dS dT , is monotone. At this point we appeal to a gener- 
alized variant of Theorem of Rado-Kneser-Choquet [15], see Theorem 2.11, 
which ensures that the harmonic extension of /: d§ dT, denoted by 
h: S T is a diffeomorphism. Certainly ^ is a quadrilateral map, i.e. 
h G Q2(S, T). We shall now see that h is an energy-minimal map. Indeed, 
by the Dirichlct energy principle and because / is a weak W^''^ -limit of fj, 
the energy of h is estimated by 

^ ^x[/] ^ liminf <fx[/j] = m. 

But m was the infimum energy in the class Q2(§, T) , so Sx[h] = m . Thus 
/iGQ2(S,T) is the energy-minimal map. □ 

The reader may see that the sequence fj actually converges strongly in 
>ri'2(S) and that f = hmS. 

Lemma 12.6. The energy-minimal map h: S T satisfies the Hopf- 
Laplace equation 

Q 

(12.1) hzhz = —2 for some c G M . 

Proof. This is immediate from Proposition 3.10, because /i is a critical point 
for all quadrilateral variations in S . □ 

Our next step is to extend ^ : S ^ T to the whole annulus X . First we 
reflect S about its horizontal side (l,i?) C 88, to obtain right-half of the 
annulus 

{z: 1 < \z\ < R, Rcz > 0} , 

and then about the imaginary axis to obtain the whole annulus X. These 
reflections, when applied to h , give rise to a mapping, denoted by ^: X 
Y , which satisfies the Hopf-Laplace equation 

(12.2) hzhz = 4 in X 

Recall that h : S T was obtained as a limit of homeomorphisms fj : S 
T converging uniformly and strongly in #^^'^(8, T) . With the aid of the 
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same reflections we obtain a sequence of liomeomorphisms fj : X Y con- 
verging uniformly and strongly in >r^'2(X,Y) to n. Thus he (X,Y). 
Now Theorem 1.4 tells us that h is an energy-minimal map within the class 
Jf2(X,Y) , and it is unique up to the rotation of X , by Theorem 1.3. In 
particular h is unique once we normalize it by setting ?i(l) = 1 . In ad- 
dition, the uniqueness combined with the fact that the quadrilateral map 
H: S T takes corners into corners in the respective order imply the 
equation 

(12.3) h{iz) = in{z) , for ah z € X. 

Now let us look more closely at the case R° , with R° small enough so 
the cracks in X are present. This yields, via Theorem 7.1, that the constant 
c in (12.4) is negative; that is, the differential hzHz dz(d)dz is positive along 
the boundary circles. According to Lemma 7.6 the cracks must lie in the 
rays of the annulus X . Since H : E> T is a diffeomorphism the cracks can 
occur only along the rays I, il, —I, — il , where I is a straight line segment 
I = [1, o] C [1, -R] ; they cannot reach the outer boundary of X because Y 
is convex at the outer boundary, see Theorem 1.6 . Moreover, in view of the 
equation (12.3), all the cracks have the same length. The proof of Theorem 
12.2 is complete. 

□ 

12.2. Lipschitz regularity, Proof of Theorem 12.3. 

12.2.1. Lipschitz regularuity of h. Local Lipschitz continuity of the admis- 
sible solutions to the Hopf-Laplace equation (with analytic right hand side) 
has been shown in [14]. Thus, it only remains to establish Lipschitz con- 
tinuity of h near the boundary points. We follow the classical method of 
extending h beyond the boundary of X in order to take advantage of the 
interior regularity results. It is characteristic of this method that sometimes 
the delicate structure of the equation is lost beyond the boundary. And that 
is exactly what happens here. For this reason we shall need a generalization 
of [14] that has been settled in [25] as follows: 

Theorem 12.7. Let h G #^^'^($7) be a mapping with nonnegative Jacobian. 
Suppose that the Hopf product hz hz is bounded and Holder continuous. 
Then h is locally Lipschitz (but not necessarily 'lo^- smooth). 

Recall that h: X Y is continuous and |^(x)| = p for jxj = R. We 
extend h to the annulus A = A(l,i?^) by the rule 

h{z) = p"^ /h{R^/z), for R^\z\<R^ 

This extension certainly belongs to the Sobolev space W^'^{A) and it has 
nonnegative Jacobian determinant. Its Hopf product equals 
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(12.4) h.h. 



' cz-"^ for Klzl^i? 



cp^ \h{F?lz)\-^z-'^ for i? ^ |z| < i?2 



Prom here we grasp the pointwise inequality < |c| . Since the Jaco- 

bian is nonnegativc, it follows that ^ l^z^z] ^ |c| , so G ^°°(A). 

Hence fiz G -^xoc i^) ' every exponent 2 ^ p < oo . Actually, hz G 
BMOioc(A) , which places H somewhere close to Lip^^ (A) but not quite 
there. Nonetheless, Sobolev imbedding theorem tells us the h is locally 
Holder continuous of exponent q; = 1 — | >0,sois the Hopf product 

hzhz ■ By Theorem 12.7 we conclude that h € Lip^^^{A) , as desired. 
Similarly, one obtains Lipschitz continuity of H near any of the four cir- 
cular arcs in dXj , which H takes into a straight line segment of dYj , by 
analogous reflection trick. 

But dealing near the corners 1, i, —1, —i G dXj requires an additional 
trick. Recall the length a of the crack, 1 < a < R and consider the sector 

CI = {z; ^ < \z\ < a , Rez > 0} = S!^(f/2 (a' ^ ^ neighborhood of 
the point 1 G dXj . We shall extend h onto CI using two different reflection 
procedures. 

• reflection in the X -space, along the arc {z = e*^ ; < ^ < 7r/2 } C 
dX , followed by the reflection in the Y -space, along the straight 
line that passes through the segment (—1, i) C dY . 

• reflection in the X -space, along the arc {z = e*^ ; — 7r/2 < 9 < 
} C 9X , followed by the reflection in the Y -space, along the 
straight line that passes through the segment {—i, 1) C dY . 

Precisely, we set the values of ^ for z = te*^ , with — 7r/2 < 9 < tt/2 and 
^ < i ^ 1 , by the rule 



(12.5) hzh-z 



-iH{l/z) + l+i , for ^ 6* < 7r/2 
h{z) = 1 , for 9 = 



in{l/z) + l-i ,for0^^<7r/2 

This is certainly a continuous map of Sobolev class W^'^{Cl,) , with non- 
negative Jacobian, and its Hopf product equals 

Q 

fizfiz = ^ ) almost everywhere in CI . 
z^ 

By Theorem 12.7 we conclude that h, G Lzpioc($^) , as desired. 

13. Two EXAMPLES 

Let us illustrate how Theorem 1.4 works in practice for mappings h : X 
Y between doubly connected domains. Thus we shall look at the Hopf dif- 
ferential \:)zi)z dz (g) dz to check as to whether it is real along dX . The two 
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examples here also serve to show a delicate difference between Hopf differ- 
entials being positive or negative. We can, and do, assume without affecting 
the results that the domain X is an annulus . Thus in either case we are 
dealing with admissible solutions to the Hopf-Laplace equation 

Q 

(13.1) ijz'^z = ^ , in an annulus X 

13.1. Case c > 0, Hopf differentials are negative along dX. This is 
the case in which no cracks emerge, by virtue of Theorem 7.1. Consider the 
following infinite series of orthogonal harmonic functions 

2 , , , R"^ - 1 ^ z"" - z-"" 

^[Zj = -T^ log U + ;r > 

(1X2) " h 

The series converges in the closed annulus A == {z; ^ \z\ ^ R} , 
except for two boundary points z = R^^. Her the symbol Log stands 

for the continuous branch of logarithm in C+ = G C ; Re ^ > } 
that is specified by Logl = 0. Observe that the expression ^ = ^^S^ 
takes values in C-|- , whenever R^^ ^ |z| ^ i? and z ^ R^^ . Clearly 
t){z) =0 for |z| = 1 and we have the following identity f)(l/z) = — f}(z) . 
Elementary geometric considerations show that t) takes the open annulus 
Ar = {z -,1 < \z\ < R} homeomorphically into a simply connected domain 
with puncture at the origin which is contained in a horizontal strip 

(13.3) f)(Afl)c{C; |ImC|<|(i?--^) } 

Passing to the limit as — )■ oo , we obtain harmonic homeomorphism 1)°° = 
z — z~^ of Aoo = ; 1 < l-zl < oo } onto the punctured complex plane 
Co . Next observe that the function ^ {z) = ^^S^ agrees with the Mobius 
transformation ^(z) = ^^S^ when restricted to any circle Cp = {z ; \z\ = 
p}. Thus the image of Cp under ^ = ^(z) is a circle in C-|- . We now 
observe that Log takes circles in C-|- into strictly convex smooth Jordan 
curves. The curves f)(Cp),l < p < R , resemble a family of ellipses with 
common focus at the origin. But the image of the outer circle, p = R, looks 
more like a parabola, but it has been flattened to fit into the horizontal strip 
(13.3). 



13.1.1. Hopf-Laplace Equation. The Hopf differential i)zi)z dz (8) dz is real 
and negative on every circle Cp , R^^ < /5 < M. These circles are horizontal 
trajectories while rays are vertical trajectories. Precisely, we have 

(13.4) i)Ji^ = X , for ah R'^ ^\z\^R, except for z = 
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Figure 3. The map i) is energy-minimal because l)zhz 

Indeed, a straightforward differentiation shows that 

Rz-1 1 , ^ R- z 1 



(13.5) 



and hi 



R — z z ^ ^ Rz — 1 z ' 

whence the equation (13.4). The Jacobian determinant 

{R' - l){\z\' - I) 



J(z,f)) = |f),|2-|f,,-|2 



\R- zl^lRz - IP 



changes sign when crossing the unit circle. 



13.1.2. The map f) is energy-minimal. In general, finding the energy-minimal 
homeomorphism between designated domains is not a trivial matter. Some- 
times it comes unplanned, like in the above example, in which the Hopf 
equation (13.4) combined with Theorem 1.4 yields: 

Proposition 13.1. Denote by Ap = f)(Ap) , where 1 < p < R . Among all 
homeomorphisms f : Ap Ap the minimum Dirichlet energy is attained 
for / = f), uniquely up to a rotation of Ap. The minimum of energy equals 



/, ,9 , ,9\ 47r log p 2tt(R^ 

{\hz\^ + lb I') = ^ + ^ 



i?2 



R^ 



1)2 . i?2 2 _ ^ 

— log ' 



R^p^ - / 



Proof. That f) represents the energy-minimal deformation and is unique up 
to the rotation follows from Theorems 1.4 and 1.3. As for the computation 
of the minimum energy, we write 



(13.6) 



i<kl<p 



Rz - 1 
R-z 



+ 



R-z 
Rz - 1 



dz 
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We expand the holomorphic functions and into Laurent series 

and use orthogonahty of the system {2;"'}n=o,±i,±2,... to arrive at the claimed 
formula. □ 

Remark 13.2. The mapping f) also represents the energy- minimal deforma- 
tion of any sub-annulus A^j \ A^^ onto a doubly connected shell \ An 
where 1< n < rs < M . The energy ^a.^\a.^ [W = '^A^, M " '^Ar, M ■ 

13.1.3. A dual Example. Somewhat dual harmonic map can be assembled 
from the orthogonal basis of harmonic functions in the annulus {z;l ^ \z\ ^ 
R } with a cut along the interval [0, R] . View it as degenerate quadrilateral. 



r r / X ^« A R^ -lyr^ 

i) = i)R{z) = -— Aigz + —— 2^ 



R ^ nR"^ 

(13.7) n=l 

2i ^ i?2-l (R-z)(zR-l) ^ , 
= -— Argz J^^^^ Jl^2 0<Argz<27r 

Straightforward computation reveals that 
We observe that the Hopf differential 

dz (8) dz 



y2 



(13.8) f)^t)^dz®dz = 

is positive along the boundary circles and negative along the cut; this is the 
energy-minimal quadrilateral deformation. 

13.2. Case c < , Hopf differentials are positive along 8%. . This 

is the case in which cracks may, though need not, emerge. We shall ex- 
ploit the typical example from the studies of the energy minimization prob- 
lem via the concept of free Lagrangians [27]. As suggested by their name, 
free Lagrangians are nonlinear differential expressions defined for mappings 
h : X Y whose integral means are independent of the boundary values, 
they depend only on the homotopy class of the mapping. These are special 
null Lagrangians. The latter were proposed by J. Ball [6] and successfully 
exploited in [10] . The most obvious novelty of free Lagrangians is that, once 
they are identified properly for a specific energy-minimization problem, the 
problem is solved instantly. 

Theorem 1.4 offers another approach. The utility of this approach is illus- 
trated by the quick proof of the following 



Theorem 13.3. Let X = A(r, R) and Y = A(r*, i?*) be planar annuli. 
Case 1. If 
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then the harmonic homeomorphism 

f,(z) = ^(- + "), {,:X°-^Y 
2 \r zJ 

has the smallest energy among all homeomorphisms h : X Y , and is 
unique up to a rotation of A. 
Case 2. If 

(13.0) -<5(-+; 

2 \ r R 

then the infimum energy among all homeomorphisms h: X '^^^ Y is not 
attained. Let a radius r < a < R be determined by the equation 

~ 2\a^ R 



r*pT r <\z\^ (J cracks along the rays [r,p\ e*^ 



Then the following mapping 

z 

^(^ + |) a ^ \z\ < R harmonic diffeomorphism 

has smallest energy within the class ,J^(K,Y) . This energy-minimal map is 
unique up to a rotation of X . 

Proof. The proof is immediate from Theorem 1.4 once we notice that 

(13.11) ^z^z dz 1^ dz = f K — , in either case 

4 

□ 
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13.2.1. The upper halves of the annuli. Let X+ = XnC"*" and Y+ = YnC"'" 
denote the upper halves of the annuli. One can view them as quadrilater- 
als with corners at {—R, —r, r, R} and {— i?* , — , , i?*} , respectively. 
Every quadrilateral map h : Y"*" extends, homeomorphically by 

reflection, into H : X Y . The energy computation, 

shows that the same mapping : X+ Y+ is the unique energy-minimal 
among strong limits of all quadrilateral homeomorphisms h : X"*" . 
Thus, in Case 2, the same cracks develop along the rays [r, p] e* ^ , < ^ < vr . 

Remark 13.4. A priori these results do not rule out the existence of univalent 
harmonic mappings from X onto Y , or quadrilateral harmonic mappings 
from X"*" onto Y+ , simply because harmonic homeomorphisms need not 
minimize the energy. Such mappings are admissible solutions to the Hopf- 
laplace equation, good enough to include them into a study of elastic defor- 
mations. Nevertheless, nonexistence of harmonic homeomorphisms between 
annuli X and Y in Case 2 was conjectured by J. C. C. Nitsche [45], [46] , 
prominent conjecture indeed. After several efforts in [17, 33, 39, 38, 40, 51], 
the conjecture was finally confirmed in [19]. Similar cracking phenomena 
are observed in higher dimensions [27], see [31] for the =Sf^ -setting. 

We end this section by showing, via a normal family argument, that 

Proposition 13.5. There is no quadrilateral harmonic homeomorphism 

h : X^ Y"*" between upper halves of annuli X = A(r, R) and Y = 

A(r*,iZ*) if the ratio ^ is too small relative to ^ . 

Proof. Fix i? > r = = 1 while varying R^ = pk > 1 , pk 1 . Suppose 
that, contrary to our claim, there exist harmonic quadrilateral mappings 
hk = Uk + ivk : X+ . The key fact is that the imaginary parts = 

Vkix, y) extend continuously up to the horizontal sides of the quadrilateral. 
Indeed, 

\m\Vk{x,£) = lim dist[/ife(a;, s), 9Y+ nM] = 

Thus the odd extension Vk{x, —y) = —Vk{x,y) beyond the horizontal 
sides defines a harmonic function, again denoted by Vk{x,y) , in the entire 
annulus X . The real part, however, may not admit harmonic extension. 
Note that each of the two sequences of harmonic functions hk : X+ Y^ 
and Vk ■ X [—Pk^ Pk] ? being bounded, contains a subsequence con- 

verging uniformly on compact subsets of X'^ and X , respectively. Let 
h: X"*" — )• C and u : X — )• M be their limits. Certainly, v is not con- 
stant in X+ , because it assumes all values from the interval (0, 1) . On the 
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other hand = 1 on X . Hence 

Q'i d — 

= TT — ^ \h{z)?' = TT (hzh + hhz) = l/izP + l/izP , so h = const, 
oz oz oz ^ ' 

clear contradiction. □ 

However, the question as to whether the same condition (13.9) draws 
the line between the existence and nonexistence of quadrilateral harmonic 
mappings h : y+ remains unclear. 

13.3. Open Question. 

Question 13.6. It is not difficult to see that if the class of harmonic diffeo- 
morphisms h : X Y is not empty, then it contains the one with smallest 
energy. Does such harmonic diffeomorphism represent the energy-minimal 
mappings among all homeomorphisms / : X Y ? The same question is 
pertinent to p -harmonic homeomorphisms. 
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